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Current vertex correction (CVC), the back-flow-like correction to the current, comes from conser¬ 
vation laws, and the CVC due to electron correlation contains information about many-body effects. 
However, it has been little understood how the CVC due to electron correlation affects the charge 
transports of a correlated multiorbital system. To improve this situation, I studied the inplane 
resistivity, p a b, and the Hall coefficient in the weak-field limit, Ah, in addition to the magnetic 
properties and the electronic structure, for a t2 9 -orbital Hubbard model on a square lattice in a 
paramagnetic state away from or near an antiferromagnetic (AF) quantum-critical point (QCP) in 
the fluctuation-exchange (FLEX) approximation with the CVCs arising from the self-energy (E), the 
Maki-Thompson (MT) irreducible four-point vertex function, and the main terms of the Aslamasov- 
Larkin (AL) one. Then, I found three main results about the CVCs. First, the main terms of the 
AL CVC does not qualitatively change the results obtained in the FLEX approximation with the 
E CVC and the MT CVC. Second, p a b and Ah near the AF QCP have high-temperature region, 
governed mainly by the E CVC, and low-temperature region, governed mainly by the E CVC and 
the MT CVC. Third, in case away from the AF QCP, the MT CVC leads to a considerable effect 
on only Ah at low temperatures, although Ah at high temperatures and p a b at all temperatures 
considered are sufficiently described by including only the E CVC. Those findings reveal several 
aspects of many-body effects on the charge transports of a correlated multiorbital system. I also 
achieved the qualitative agreement with several experiments of Sr2RuC>4 or Sr2Ruo.975^0.02504. 
Moreover, I showed several better points of this theory than other theories. 

PACS numbers: 71.27.+a, 74.70.Pq 


I. INTRODUCTION 

Many-body effects, effects of Coulomb interaction be¬ 
tween itinerant electrons beyond the mean-field approxi¬ 
mation, are important to discuss electronic propertied 14 . 
When the Coulomb interaction is very small compared 
with the bandwidth of itinerant electrons, which is of 
the order of magnitude 1 eV, we can sufficiently de¬ 
scribe its effects in the mean-field approximation as the 
static and effectively single-body potentiald^l However, 
in correlated electron systems such as transition met¬ 
als or transition-metal oxides, the Coulomb interaction 
becomes moderately strong or strong, resulting in the 
derivations of their electronic properties from the single¬ 
body picture 1 ®. 

For several correlated electron systems, many-body ef¬ 
fects can be described in Landau’s Fermi-liquid (FL) the¬ 
ory' 10 . This theory is based on two basic assumptions?®. 
One is the one-to-one correspondence between the non¬ 
interacting and the interacting systems. Because of this 
assumption, we can describe low-energy excitations of the 
interacting system in terms of quasiparticles (QPs) with 
the renormalized effective mass and the renormalized in¬ 
teractions described by the Landau parameters 10 . The 
other assumption is lack of the temperature dependence 
of the Landau parameters. Because of this assumption, 
the temperature dependence of the electronic properties 
remains the same as that in the noninteracting systenP®. 
Furthermore, as a result of those assumptions, many- 
body effects on the electronic properties are the changes 
of their coefficients due to the mass enhancement or the 


FL correction or botiM 

Actually, Landau’s FL theory well describes several 
electronic properties of S^RuOq at low temperatures. 
First, this theory can explain the almost temperature- 
independent spin susceptibility 11 and the T 2 dependence 
of the inplane resistivity 12 . In addition, the importance 
of many-body effects has been suggested in the measure¬ 
ments of the de Haas-van Alphen (dHvA) effect 18 ’ 1 ^ and 
the Wilson raticP-? the effective mass of the d xz / yz or 
the d xy orbital measured in the dHvA becomes, respec¬ 
tively, 3 — 3.5 or 5.5 times as large as the mass obtained 
in the local-density approximation (LDA), a mean-field- 
type approximation; the Wilson ratio, the ratio of the 
spin susceptibility to the coefficient of the electronic spe¬ 
cific heat, becomes 1 . 7 — 1.9 times as large as the noninter¬ 
acting value. Note that the enhancement of the Wilson 
ratio arises from the FL correction*^. 

However, we observe non-FL-like behaviors, the de¬ 
viations from the temperature dependence expected in 
Landau’s FL theory, for correlated electron systems near 
a magnetic quantum-critical point (QCP)P®. For exam¬ 
ple, Sr2Ruo.975Tio.o25 04, a paramagnetic (PM) ruthen- 
ate near an antiferromagnetic (AF) QCP, shows the 
Curie-Weiss-type temperature dependence of the sp in 
susceptibility and the T-linear inplane resistivity * 16 * 17 ! 
Also, Ca2- x Sr r Ru04 around x = 0 . 5 , a PM ruthen- 
ate near a ferro magne tic QCP, shows the similar non- 
FL-like behavior d 18 * 19 *. Thus, those experimental results 
indicate the importance of many-body effects beyond 
Landau’s FL theory near a magnetic QCP. Note, first, 
that the wave vector of the spin susceptibility enhanced 
most strongly in S^Ruq. 975^0.025O4 20 ‘ is the same for 
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Sr 2 RuOi 2 H, i.e. q « (^, ^); second, that Ti substi¬ 
tution does not cause any RuOg distortion^®, while Ca 
substitution causes RuOg distortions such as the rotation 
and the tiltin pP^, which drastically affect the electronic 
structur e) 23 ! 24 ). 

Among correlated electron systems, the ruthenates are 
suitable to deduce general or characteristic aspects of 
many-body effects in correlated multiorbital systems be¬ 
cause of the following three advantages. The first advan¬ 
tage is that the ruthenates show the FL or the non-FL- 
like behaviors, depending on the chemical composition 
or the crystal structure or botlf 1112 fSHH] Due to this 
advantage, we can study how the FL state is realized 
and how the system changes from the FL state to the 
non-FL-like state, and we may obtain their general or 
characteristic properties. Then, the second advantage is 
that the ruthenates are the t 2o -orbi tal sy stems with mod¬ 
erately strong electron correlatior l 14 l 25 [ This has been 
established for Sr 2 Ru 04 by three facts: the Ru t 2g or¬ 
bitals are the main components of t he den sity-of-s tates 
(DOS) near the Fermi level in the LDA^ 2 ^ the LDA^ 27 ) 
can reproduce the topology of the Fermi surface (FS) ob¬ 
served experimentally 1 - 3 23 the experimentally estimated 
value of U , onsite intraorbital Coulomb interaction, is 
about 2 eV^l, which is half of the bandwidth for the 
t 2g orbitals in the LDAEi3 2 3. In addition to the second 
advantage, the third advantage is the simple electronic 
structures ' compared with the other multiorbital sys¬ 
tems. Due to the second and the third advantage, we can 
simply analyze many-body effects of a correlated multi¬ 
orbital system, and that analysis may lead to a deep un¬ 
derstanding of the general or characteristic aspects of the 
many-body effects. 

To describe many-body effects near a magnetic QCP, 
we need to use the theories that can satisfactorily take 
account of the effects of the critical electron-hole scat¬ 
tering arising from the characteristic spin fluctuation of 
that QCP. If the system approaches a magnetic QCP, we 
observe the enhancement of the spin fl uctuat ion for the 
wave vector characteristic of that QCFL1M29I That en¬ 
hancement causes the strong temperature dependent crit¬ 
ical electron-hole scattering mediated by the spin fluctu¬ 
ation. Then, that critical electron-hole scattering results 
in the emergence of both the hot spot of the QP damp¬ 
ing and the Curie-Weiss-type temperature dependence of 
the reducible four-point vertex function for the momenta 
connected by the spin fluctuation. (Note that the re¬ 
ducible four-point vertex function describes the multiple 
electron-hole scatterings).) Thus, the emergence of the 
former violates the first basic assumption of Landau’s 
FL theory since at the hot spot the QP lifetime is not 
so sufficiently long as to realize an approximate eigen¬ 
state as Landau’s FLP. Furthermore, the latter violates 
the second basic assumption because the temperature de¬ 
pendence of the reducible four-point vertex function and 
mass enhancement factor determines the temperature de¬ 
pendence of the Landau parameter) 10 . Thus, many-body 
effects near a magnetic QCP may be described by the 


theories beyond Landau’s FL theory if the theories can 
satisfactorily treat the strongly enhanced temperature- 
dependent spin fluctuation. 

Actually, several non-FL-like behaviors near a mag¬ 
netic QCP can be reproduced in fluctuation-exchange 
(FLEX) approximatiorPSED with the current vertex cor¬ 
rections (CVCs) arising from the self-energy (X) and 
the Ma ki-Th ompson (MT) irreducible fo ur-po int vertex 
functioii 35 ) 36 ! due to electron correlatior l 37 ) 38 ) For ex¬ 
ample, this theory shows the Curie-Weiss-type tempera¬ 
ture dependence of both the spin susceptibility and the 
Hall coefficient and the T-linear inplane resistivity for 
a single-orbital Hubbard model on a square lattice in a 
PM state near an AF QCP, where the spin fluctuation 
for q = ( 7 r, 7 r) is enhanced 37 . Those results are consis¬ 
tent with the experiments of cuprates®)®). Since the 
powerfulness of the FLEX approximation near a mag¬ 
netic QCP arises from its satisfactory treatment of the 
momentum and temperature dependence of spin fluctu- 
ationfhHHD the similar applicability will hold even for a 
multiorbital Hubbard model on a square lattice. 

Since the effects of the CVCs due to electron-electron 
interaction in a correlated multiorbital system had been 
unclear, I studied several electronic properties of an ef¬ 
fective model of several ruthenates, a t 2s -orbital Hub¬ 
bard model on a square lattice in a PM state away 
from or near an AF QCP, in the FLEX approximation 
with the X CVC and the MT CVCPEHl, and then I ob¬ 
tained satisfactory agreement with several experiments 
and three important aspects of many-body effects on the 
charge transports. First, the results away from the AF 
QCP qualitatively agree with five experimental results 
of Sr 2 Ru 04 , (i) the strongest enhancement of spin fluc¬ 
tuation 21 at q ss (^, ^), (ii) the nearly temperature- 
independe nt sp in susceptibility®, (iii) the larger mass en¬ 
hancement - 3 14! -of the d xy orbital than that of the d xz / yz 
orbital, (iv) the T 2 dependence of the inplane resistiv¬ 
ity at low temperatures®, and (v) the non-monotonic 
temperature dependence of the Hall coefficient®. Note 
that the Hall coefficient observed in Sr 2 RuC >4 shows the 
following non-monotonic temperature dependenc^E a t 
high temperatures above 130K, the Hall coefficient is 
small and negative with a slight increase with decreas¬ 
ing temperature; after crossing over zero at 130K, the 
Hall coefficient becomes positive with keeping an in¬ 
crease, and shows a peak at about 70K; below about 
70K, the Hall coefficient monotonically decreases with 
decreasing temperature. Then, the results near the AF 
QCP can qualitatively explain three experimental results 
of Sr 2 Ruo. 975 Tio.o 25 C> 4 , (i) the strongest enhancement of 
spin fluctuation!^ at q ~ (yy,yy), (ii) the Curie-Weiss- 
type temperature dependence of the spin susceptibility®, 
and (iii) the T-linear inplane resistivity 17 . In this com¬ 
parison, I assume that the main effect of Ti substitution 
is approaching the AF QCP compared with Sr 2 Ru 04 . 
Note that the measurement of the Hall coefficient in 
Sr 2 Ru 0 . 9 _ 75 Ti 0 . 02 . 5 O 4 has been restricted at very low tem¬ 
perature 43 ), which is out of the region I considered. More- 
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over, I revealed the realization of the orbital-dependent 
transports, the emergence of a peak of the temperature 
dependence of the Hall coefficient, and the absence of the 
Curie-Weiss-type temperature dependence of the Hall co¬ 
efficient near the AF QCP. 

However, the previous studied 121221 contain two remain¬ 
ing issues. One is to clarify many-body effects of the 
Aslamasov-Larkin (AL) CVC, the CVC ar ising from the 
AL irreducible f our-p oint vertex functiorPESl. l n the 
previous studies 1 ®^ i neglected the AL CVC in the 
FLEX approximation for simplicity since in a single¬ 
orbital Hubbard modeP^on a square lattice the AL CVC 
does not qualitatively change the results of the resistivity 
and Hall coefficient near an AF QCP and since the sim¬ 
ilar property would hold even in a multiorbital Hubbard 
model on a square lattice not far away from an AF QCP. 
However, it is necessary and important to analyze the ef¬ 
fects of the AL CVC in that multiorbital Hubbard model 
since both the MT and the AL CVC are essential to hold 
conservation laws exactly 1 ®. In particular, that analysis 
is needed not only to check the validity neglecting the AL 
CVC for qualitative discussions but also to clarify many- 
body effects of the AL CVC. The other remaining issue is 
to give the comprehensive explanations about the formal 
derivations both of the transport coefficients in the ex¬ 
tended Eliashberg theory®! to a multiorbital system and 
of the CVCs in the FLEX approximation for a multior¬ 
bital Hubbard model. My previous study^l reported a 
microscopic study about the effects of the CVCs due to 
electron correlation in a multiorbital system. In the pre¬ 
vious studiephEH], however, I just gave brief explanations 
about those formal derivations. Thus, it is desirable to 
explain the detail of those formal derivations since those 
will be useful to adopt the same or similar method to the 
transport properties of other correlated electron systems. 

In this paper, after formulating the dc longitudinal 
and the dc transverse conductivities in the extended 
Eliashberg theory to a multiorbital Hubbard model in the 
FLEX approximation with the CVCs, I study the effects 
of the main terms of the AL CVC on the in-plane resistiv¬ 
ity and the Hall coefficient for the quasi-two-dimensional 
PM ruthenates near and away from the AF QCP. As the 
main results, I show the quali tative validity of the main 
results of the previous studies 1 ®^, the existence of two 
almost distinct regions of the charge transports near the 
AF QCP as a function of temperature, and the different 
effects of the MT CVC on the low-temperature values of 
the in-plane resistivity and Hall coefficient away from the 
AF QCP in the presence of the X CVC, the MT CVC, 
and the main terms of the AL CVC. I also present sev¬ 
eral results of the magnetic properties and the electronic 
structure, and show four main results about each of the 
magnetic properties and the electronic structure. Those 
are useful for deeper understanding than in the previous 
studied®. 

The remaining part of this paper is organized as fol¬ 
lows. In Sec. II, I explain the method to calculate the 
electronic properties of some quasi-two-dimensional PM 


ruthenates without the RuC>6 distortions. In Sec. II A, 

I show the Hamiltonian of an effective model of some 
quasi-two-dimensional ruthenates, explain the parameter 
choice for the noninteracting Hamiltonian, and briefly re¬ 
mark on the spin-orbit interaction. In Secs. II B 1 and 

II B 2 , I explain the extended Eliashberg theory to the 
dc longitudinal and the dc transverse conductivities for a 
multiorbital system and give several theoretical remarks 
about their general properties. In Sec. II C, I explain 
several advantages of the FLEX approximation with the 
CVCs, formulate the FLEX approximation for a multi¬ 
orbital Hubbard model, and derive the Bethe-Salpeter 
equation for the current with the X CVC, the MT CVC, 
and the AL CVC in the FLEX approximation. Further¬ 
more, I derive a simplified Bethe-Salpeter equation by 
approximating the AL CVC to its main terms. In Sec. 
Ill, I show the results of several electronic properties of 
the quasi-two-dimensional PM ruthenates near and away 
from the AF QCP in the FLEX approximation with the X 
CVC, the MT CVC, and the main terms of the AL CVC; 
in addition to that case, I consider three other cases con¬ 
sidered in Ref. m for discussions about the transport 
properties in order to deduce the main effects of the AL 
CVC. After discussing the magnetic properties in Sec. Ill 
A and the electronic structure in Sec. Ill B, I discuss the 
main effects of the AL CVC on the inplane resistivity and 
the Hall coefficient in Sec. Ill C. Then, I compare the 
obtained results with several experiments of S^RuCH or 
Sr2Ru0.975Ti0.02.5O4 in Sec. IV A, and other theories in 
Sec. IV B. In Sec. V, I summarize the obtained results 
and their conclusions, and show several remaining issues. 


II. METHOD 

In this section, I explain an effective model of some 
quasi-two-dimensional ruthenates and a general theoret¬ 
ical method to analyze the resistivity and the Hall coeffi¬ 
cient for a correlated multiorbital system in a PM state. 
In Sect. II A, we see the Hamiltonian of the effective 
model, determine the parameters of the noninteracting 
Hamiltonian, and remark on the spin-orbit interaction, 
neglected in the effective model. In Sect. II B 1 , to ana¬ 
lyze the resistivity, we explain the formal derivation of the 
dc longitudinal conductivity of a multiorbital Hubbard 
model in a PM state without an external magnetic field 
in the linear-response theory with the most-divergent- 
term approximation^, and show general properties of 
the derived longitudinal conductivity and their conse¬ 
quences for the resistivity. In Sect. II B 2 , we derive 
the dc transverse conductivity of a multiorbital Hubbard 
model in a PM state in the weak-field limit by using the 
linear-response theory with the most-divergent-term ap¬ 
proximation, and see general properties of the derived 
transverse conductivity and the Hall coefficient in com¬ 
bination with the results for the longitudinal conductiv¬ 
ity. The general formulations in Sects. II B 1 and II 
B 2 are the extensions of the single-orbital cases for the 
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FIG. 1: Schematic pictures of the hopping processes of the 
t 2 g orbitals on a two-dimensional layer. The difference in the 
color of each orbital shows the difference in the sign of its 
wave function. 


resistivity"® and the Hall coefficient 4 !!®, respectively. In 
Sect. II C, we remark on several advantages of the FLEX 
approximation with the CVCs, formulate the FLEX ap¬ 
proximation in Matsubara-frequency representation for a 
multiorbital Hubbard model in a PM state in the similar 
way for Refs. :33] and SH and derive the CVCs in the 
FLEX approximation by extending the formulation for a 
single-orbital case®. 

Hereafter, we use the following unit and notations: We 
set 7i=l,c=l,e = l, /xb = 1, and k-Q = 1. In the equa¬ 
tions, the d xz , d yz , and d xy orbitals are labeled 1 , 2 , and 
3, respectively. In Matsubara-frequency representation of 
several quantities, we use the fermionic and the bosonic 
Matsubara frequency, e m = irT(2m+l) and f l n = 2 ttT n, 
respectively. In real-frequency representation, we use 
frequency variables such as e and to and abbreviations 
such as k = (fc, e) and q = ( q,u )), with momenta k and 
q. We use the abbreviations such as X{a} — b c tP 
S{si} = XX,s2,s 3 ,s 4 ’ an d r { a }{k,ie m ,k , ie m '] 0, i£l n ) = 
Tabcd(fc, id m-> k i 


A. Effective model 

In this section, I introduce the total Hamiltonian of an 
effective model for some quasi-two-dimensional ruthen- 
ates and explain how to choose the parameters in the 
noninteracting Hamiltonian. I also give a brief remark 
about the spin-orbit interaction. 

To describe the electronic properties of several 214- 
type rutlre nates such as S^RuCH, I use a ^g-orbital Hub¬ 
bard modet®® on a square lattice because several 214- 
type ruthenates are categorized as quasi-two-dimensional 
^ 2 g _ orbital correlated systems and Ru ions on a two- 
dimensional layer form a square lattice 14 . The Hamil¬ 
tonian of this model is 

H = H 0 + Hint, (1) 

where Hq and H; nt are the noninteracting and the inter¬ 
acting Hamiltonian, respectively. 

First, Ho is given by 

3 

Hq = e a6(fc)4as^bs- (^) 

k a,6—1 s=t>4 


Here Cfc as and c^. as are the annihilation and the creation 
operator of an electron of momentum fc, orbital a, and 
spin s, and e a b(k) is given by 


bi(fc) — 

A *2 S 

3 

— 2 ti cos k x 

— 2 1 2 cos k y — /x, 

(3) 

£ 12 (fc) = 

£21 (fc) 

= At' sin k x 

sin k y , 

(4) 

£22 (k) = 

A *2 3 

3 

— 2 f 2 cos k x 

— 2f 1 cos k y — /x, 

(5) 


2A t . 




£ 33 (fc) = 

3 

— 2t 3 (cosfc x 

+ cos k y ) 



— 4^4 cos k x cos k y - 

~ /b 

(6) 


and otherwise e a b{k) = 0 , where A t2 is the difference 
between the crystalline-electric-field energies of the d xy 
and the d xz / yz orbital, p, is the chemical potential deter¬ 
mined so that the electron number per site, n e , satisfies 
n e = 4, and t\ , t 2 , £ 3 , Li, and t' are the hopping integrals 
of the tjgn orbitals, whose schemat ic pic tures are shown 
in Fig. [4 Since I neglect the effeetd 23 ^® of the RuC >6 dis¬ 
tortions on H 0 , the targets of this paper are the 214-type 
ruthenates without the RuC >6 distortions. 

Assuming that the l.l )A— 2j for Si^RuCH gives a good 
starting point to include many-body effects in the 214- 
type ruthenates without the Ru0 6 distortions, I choose 
the parameters in e a b(k) so as to reproduce the electronic 
structure obtained in the LDA^ for S^RuCH. Namely, 
I set t\ = 0.675 eV, 1 2 = 0.09 eV, to = 0.45 eV, t± = 0.18 
eV, t' = 0.03 eV, and A t2g = 0.13 eV. As explained in 
Ref. [TOJ the obtai ned electronic structure is consistent 
with the LDA^X the bandwidth for the t 2g orbitals is 
about 4 eV; the quasi-one-dimensional d xz and d yz or¬ 
bitals form the hole-like a and electron-like /3 sheets, and 
the quasi-two-dimensional d xy orbital forms the electron¬ 
like 7 sheet; the van Hove singularity of the d xy orbital 
exists above the Fermi level; the occupation numbers of 
the d xz / yz and the d xy orbital are n xz / yz — 1.38 and 
‘blxy = 1.25. 

Then, H lnt is given by 


TT _1 rrSlS2S 3 S 4 -t "t 

Hint -A 2^ 2^ 2^ U ab, 


abed 'jas , 'jd .s4 f j rs 3 £'J bs 2 

3 {“} {si} 

3 3 

= ^££ Tlj a'f Tlj a], 4 “ U £££ Tlja'ft'jb 
j a —1 j a= 16 <a 

3 1 

— Ju + 2 ^ja^jb) 

j a— 1 b<a 
3 

+ J, £££ 

j a— 1 b^a 


(7) 


Here H*^J S3S4 is a bare four-point vertex function, U 
is intraobital Coulomb interaction, U' is interorbital 
Coulomb interaction, Jh is Hund’s rule coupling, J' is 
pair hopping term, h ja is n ja = X) s n jas = X s ^ jas c jas , 
and s ja is s ja = \ YL s ,s' ^j as a ss'C-jas’ with the Pauli ma¬ 
trices (j ss f . Among the terms of U 3 ^ S3Si , it is sufficient 
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for a PM state to use U^ cd , C/^L, and U± cd , which are, 
respectively, = U$£ = U^J, uH cd ee = 

and u abcd = U lt!d = In addition, in the 

absence of the spin-orbit interaction, it is more useful 
to introduce bare four-point vertex functions in spin and 
charge sector, U^ bcd and U^ bcd , defined as 


77 S 1 S 2 S 3 S 4 
^ abed 


1 

2 


u: 


abcd L 


0 

SlS 2 


G 


0 

S4S3 


177 S . 

> ^abed^ 


( 8 ) 


Namely, C/f bcd is U^ bcd = U^ cd + uH cd , 
Uabcd = U^bcd ~ U^lcd = —U± bcd - 


and u abcd is 


I will explain how to treat the effects of H lnt in Sec. II 
C, and how to choose the values of U, U ', Jh, and J' in 
Sec. III. 

In the effective model, I neglect the spin-orbit inter¬ 
action of the Ru t 2 g orbitals for simplicity. This treat¬ 
ment may be sufficient to discuss the electronic properties 
analyzed in this paper since the couplin g c onstant esti¬ 
mated in local-spin-density approximation 49 for Si^RuCH 
is 0.167 eV, which is smaller than the main terms of 
Ho and H mt , and since its effects will not qualitatively 
change the results shown in Sect. III. (The main terms 
of H int are of the order of magnitude 1 eV, as described 
in Sec. III.) For several expected roles of the spin-orbit 
interaction, see the remaining issues in Sec. V. 


B. Extended Eliashberg theory to charge 
transports of a multiorbital system 


Before the formal derivations, I remark on the mean¬ 
ings of taking the ui limit and holding wrt ra ns 1 in 
these derivation^- 19 with r trans , the transport relaxation 
tim d 46 * 47 ! (of the order of magnitude the QP damping). 
First, the u> limit, i.e. lim w _>-o liniq-^O) is vital to obtain 
the observable currents since the dynamic and uniform 
field causes the observable currents; on the other hand, 
the q limit, i.e. lim q ^.o lim^-^Oi does not cause any ob¬ 
servable currents as a result of the screening induced by 
the modulations of the charge distribution 50 . Then, in 
taking lim^^o, the QP lifetime should hold wr trans -C 1 
since the inequality characterizes the relaxation process 
of transports; in WTt rans -C 1, local equilibrium is realized 
due to the rapid relaxation compared with w _1 , a typical 
time scale of the field, and then the QPs near the Fermi 
level mainly govern the electronic transports. 


1. Resistivity 


For discussions about the resistivity of a correlated 
multiorbital system, I use the Kubo formulapl for the lon¬ 
gitudinal conductivity, a vv , in the co limit and tuT trans <C 
1 , 


k£\q,u))-K£>(q,0) 


a vv = 2 lim lim 

uj —>-0 q —>0 


(R)/ 


uu 

= 2 lim *ff ) (0.*)-e , (0.Q) , 

^-►0 iuj 


where Kw\0 ,w) is determined by Kw\0 ,w) = 

K uv {iVL n —> w + *0+) with K vv {iQ, n ), being 


In this section, we derive the dc longitudinal conductiv¬ 
ity without an external magnetic field and the dc trans¬ 
verse conductivity in a weak-field limit in the linear- 
response theory with the most-divergent-term approx¬ 
imation. In Sec. II B 1, we derive the dc longitudi¬ 
nal conductivity to analyze the resistivity of a correlated 
multiorbital system. After deriving the exact expression 
in terms of the four-point vertex function or the three- 
point vector vertex function, we derive its approximate 
expression in the most-divergent-term approximation"®, 
which is appropriate for the metallic systems with long- 
lived QPs at (at least) several momenta. We also explain 
four general properties seen from the derived expression 
of the conductivity and show the properties of the re¬ 
sistivity about the dominant excitations, the dependence 
on the QP lifetime, and the main effects of the CVCs. In 
Sec. II B 2, to analyze the Hall coefficient of a correlated 
multiorbital system for a weak magnetic field, we derive 
the dc transverse conductivity in the weak-field limit. 
Due to difficulty deriving the exact expression, I derive 
only the approxima te ex pression in the most-divergent- 
term approximation^ 7 In addition, after explaining 
four general properties of the derived conductivity, we 
deduce the properties of the Hall coefficient in the weak- 
field limit about the similar things for the resistivity. 


rp—1 

= lim l / d Te *Wr(T T j (r )j (o)) 

qr^O J Q 

1 rT- 1 

= ArEE/ dTe innT {v kv ) ba {v k ' u ) cd 

V k,k' {a} 170 

x ( t t 4 &(r)c fe a(r)c^ c Cfc'd) 

J 1 _ 

= -^EEE( Wfc v)ba{Vkv) cdG aci}^ i ^m-\-n)G dbi}^ i ^m) 

k m {a} 

^2 _^ ^ ^ ^ 

_ ~jtf 2 ^ v ^ y ^^^^{ykv)ba{Vk'v)cdGaA{k’>'i'€m+n) 
k,k' m,m' {a} {A} 

X GdD^fc 1 Bb^k] 5 ^m'+n) 

x T {A} (k,ie 

mi k 1 fom' j 0 ) (Id) 


In Eq. (10), ( v kv ) ab is the group velocity, 


(t J k iy)ab — 


de a b(k) 

dk v 


( 11 ) 


and P{A} (fe, ie m , k', ie m r,0 , ifl n ) is the reducible four- 
point vertex function, which is connected with the ir- 
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reducible one through the Bethe-Salpeter equation^, 


the detail see Appendix B): 


I '{A}(k,ie mi k i ie m /; 0 , ifi n ) 

= (k, ie m , k , i€ m ' 5 0 , 

T 

+ atEEE r ABC'D’{k,ie m 1 k , *e m »;0,*f2„) 

fc" m" {A'} 

X GcA'(k 5 ^m // +n)G'B / D / ?^m") 

x r^V B , CD (fe",«e m ",fc',*e m /; 0 ,i^n)- ( 12 ) 

The irreducible four-point vertex function can be deter¬ 
mined in the way explained in Sec. II C. 

To obtain an exact expr ession of ay„, we carry out the 
analytic continuations®^ 0 f the first and second terms 
of Eq. (10 1 by using the analytic properties of the single¬ 
particle Green’s function and the four-point vertex func¬ 
tion. As we will carry out those analytic continuations 
in Appendix A, we obtain 




de 
47 vi 


is Jed, 

k,k' {a} 

tanh 2 ^Eq ) a}( fc ’ fe,;e;W ^ 

(tanh - tanh /T™ } (fc, k'; e; u) 


— tanh 


2 T 


K. 


•g ) a} (fc,fe ,; e ; w) 


(13) 


with Ki^ a y (fc, fc'; e; w), being 

E{o}( fc,fe,;e;W ) = 9 l;acdb(k-,u)S ki k' 

1 /- 00 de' x - ^ 

iv / ^ ^ 


, „ 47TJ 

°° {A}i'=l 


aABb 


(fc; w) 


X x 7 «';{A} (fc, fc'; W )gv-CcdD(k w). 

(14) 

Here Jw,{A}{k,k' ;w) is connected with the reducible 
four-point vertex function in real-frequency representa¬ 
tion and is determined by the Bethe-Salpeter equation, 

J'«' ; {A}(fc,fc';w) = J^} {A} {k,k'-,u}) 

+ E AfE E / T-Jii";ABC'D'(fc,fc ! w) 


' = 1 


fc" {A'}' 


47T* 

rl 1 


x 9l";C'A’B'D'{k i ul )Jl"l'- J A'B l CD(k jfcjw), (15) 


^(O.w) 

1 

= ~~ TV 


/ o 

-( 


k {a} 


de 
47 ri 


x 


tanh ^ 9i-acdb(k ; w)A„ ;1;cd (fc; w) 


(tanh - tanh gi-,acdb{k] w)A„ ; 2 ; C d(fc; w) 

(16) 


- tanh £ ^ 93-,acdb{k', cv)A v . 3;cd (t, w) 


Because of the difficulty solving the exact expression 
of OV;,, we use the most-divergent-term approximation, 
introduced by Eliashberg®, in order to derive an ap¬ 
proximate expression. In this approximation®, we con¬ 
sider only the most divergent terms with respect to the 
QP lifetime in 7 *(fcp)/T —► 0 with 75(fcF)> the QP 
damping for band a at Fermi momentum fc F . This ap¬ 
proximation is based on the limiting properties 8 9 of the 
pairs of two single-particle Green’s functions with exter¬ 
nal momentum and frequency, q and ui, in q 0 and 
7 *(fc F )/T —► 0. More precisely, utilizing the limiting 
properties, we can use the approximation that among 
the pairs of two single-particle Green’s functions, only 
a retarded-advanced pair gives the leading dependence 
on the QP damping and the external momentum and 
frequency. Namely, we can approximate the leading de¬ 
pendence of gi-acdb(k;ui ), g 2 - ac db(k;uj), and g 3 -,acdb(k; u) 

tcP 


g\\acdb{k‘i tu) 


52;acdfc(fc;w) 


Yli U ^)aa(Ul) a c(U h ) d p(Ul)p b 

a,(3 

_ z a (k)zp(k) _ 

[ e — £a(fc) + *0+][e — Q(fc) + T0+] 

(17) 

27 xi Y,( U ^(Ut) ac (U k U(Ut) 0b 

a, (3 

Zg(k)zp(k)6(e -&(k)) 

u - C(fc) + Cp{k) + *[7j(fc) + 7|(fc)] ’ 

(18) 


and 


9^\acdb{k\ (d) ^ 52(Uk)aa(Vl)ac(Uk)df>(Vlhb 

a,f3 

_ Z a (k)zp(k) _ 

X [e - C(fc) - *°+H e - Q(fc) - * 0 +] ’ 

(19) 


with the similar connection between J^j^\ A ,[k,k'\ui) 
and the irreducible four-point vertex function in real- 
frequency representation. 

~ /r>\ 

We also rewrite Ay,„ ( 0 , to) in a more compact form by 
using the three-point vertex function in real-frequency 
representation, A v -i- a b{k; uj) = A u -i. ab (k, e + w, fc, e) (for 


respectively. Here ££(fc) is the QP energy, z a {k) is 
the mass enhancement factor, and ( Uk) aa is the unitary 
matrix to obtain the QP dispersions. Since this treat¬ 
ment remains reasonable for 'y^ l (kp)/T < 1 , the most- 
divergent-term approximation is not only exact in the FL 
but also appropriate in the correlated metallic systems 
having the long-lived QPs at least for several momenta. 
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To derive an approximate expression of ovj, in the 
most-divergent-term approximatioiP^, we introduce two 
quantities, and A^. ab (k;ui), which are 

irreducible only about a retarded-advanced pair, and 
rewrite K^/ 1 (0. cj) by using the two quantities. First, 
we define J$\ a Ak,k']U)) and A ^.}. ab {k;co) as 




(i) 


de 


^EE / h E 


k" {A}' 




l"= 1,3 
r(l) 


x gi";CABD(k 'i^lJl'iy.ABcd!^ ’k;U)), ( 20 ) 


and 


and 

1 

N 


EE 

k {a}' 


. . (Vkv)ba t&nh Q^acdbi}^i ^0 

47TZ 2i 


*ee/“ <^32;cdCD ^ > w )d2 ;CABD(k \U>) 


k' {A} 1 ' -00 
X A 1 /; 2 ;AB(fc'; w) 

1 A 50 de 


47U liv 


1 


EE 


= -^EE 

k {o} " ' fc' {A} 

r(°) ei- , i- , i . . i.' J iJ J 


00 dd 

-rv, 47 TI 


x ^ 23 ; 6 oBa(^’ e, fc, e + u), k! , e , fc', e' + w) 
x 93-,BDCA{k', e', fc', e' + w)(dh,v)dc 

x (tanh —^ 7 - ^ tanh ^^ 2 ;acdh(fc;a;)A y; 2 ; C d(fc; w), 

(24) 


A i/;Z;ab(^> W ) — { v ki^)ab 

AEjEfJiUM 

{A}Z' = 1,3 fc' 00 


x gi’-CABD(k'\u){Vk’v)AB, 


( 21 ) 


respectively. We can also connect A ( „.}. ab (k;u)) with 
A v -i-ab{k\oj) as follows: 


A v -i-ab{k\u) = A ( °l ab (k;uj) 

{A} fc' • /_0 ° 

x g 2 -CABD(k’; U))A v -2-AB(k'; w). 


( 22 ) 


Then, substituting Eq. (22) into Eq. (16) and using two 
equalities, 


-apEE / ;pr(^)&atanh -*=eq 

;acdb {^5 

fc {a} *' _0 ° ^ 

4 eeT > ■•^12-cdCD^ik ;uj)g 2 -cABD{k ;w) 

fc' {A} d ~°° 

x A„ ;2 ;As(fc 7 ; w) 

1 w f 00 de \ 1 V V 

“ ~~ TV ^ ^ I oo 4 tt 7 LiV ^ ^ ,, ^ 

fc {a} J_ °° fc' {A} J 

X ^21;6aBA(^> e ; 6 + ^ ^ + w ) 

x 9i-,bdca( k\ e', fc', e' + w)(wfc<„)flc 
x (tanh —^r - tanh -^j g 2 -acdb(k; co)A„ ;2;cd (k; w), 

(23) 


00 dd 

47TZ 


~ /r>\ 

we can express iOi/( 0 ,cj) as two parts, the part exclud¬ 
ing a retarded-advanced pair and the other part: 

(°. w ) = - t? E E^"^ f j— 


fc {a} 


tailh ^f gi ’ acdb ( k ' W ) A l;l;cd( fc ! w ) 


47T* 

(o) 

1 ; 

( 0 ) 


- tanh ~Y^93;acdb{k-, w)A(,.'. C(i (fc; u) 


-^EE 


fc {a} 


^ A S; b j fc ’ e ’ fe ’ e + W ) 


(tanh —-tanh^) 52 ;acd&(fc;w) 


x A„ ; 2 ;c d(fc;w). 


(25) 


This expression remains exact at this stage. In Eqs. (23) 


and ( 24), we have used Eqs. (|A7[), (A9|, (All), and 
(A13) and the exchange symmetry^ of the four-point 


vertex function about its variables. 

Adopting the most-divergent term approximation to 
Eq. (251, extracting the w-linear term, and using ([9]), we 
obtain an approximate expression of oy^, 


y' d e ( ^-C( e ) 

iv 2^ 2^ / l 


fc {a}=l 1 


27r V de 


X Ay° 2 ;ba(^; 92-,acdb{k\ 0) A„ ; 2 ;c d(fc; 0 ). (26) 

Here we can regard A^.\. ba (k; 0) and A„ ; 2 ; cd(fc;0) as, re¬ 
spectively, the current including the CVC arising from 
the self-energy and the current including the CVCs aris¬ 
ing from the self-energy and the irreducible four-point 
vertex function. This is because Eq. (21) for l = 2 at 
w = 0 becomes 


A l;Lb( fc ; °) = (v k „)ab + 


dReE if(k) 

dk v 


(27) 




























as a result of a Ward identitjP, and because Eq. (22 1 for 
l = 2 at u> = 0 becomes 


A^M) = A^ ;cd 



(M) 

^ T (i) 

47T* 22 ;cdCD 


(fc,fc';0) 


x 92-,CABD(k r ; 0)A„ ;2 -AB(k'; 0), 


(28) 


as a result of the disappearance of the second term of Eq. 
(20), the higher-order terirP^ about WTt ran s than the first 


term of Eq. (20). Note that the second term of Eq. (28) 
plays a similar role for the backflow correctiorP in the FL 
theory since that term connects the currents at k and k'. 


From Eq. (26 1 , we see four general properties for the 


dc longitudinal conductivity of a correlated electron sys¬ 
tem. (The following arguments for <j xx are qualitatively 
the same even for cr yy .) First, due to the factor (— 9 ^ ) 
in Eq. (26), the main excitations arise from the QPs 


near the Fermi level. This property indicates the impor¬ 
tance of the coherent part of the single-particle Green’s 
function in discussing p a {,. Such importance holds even 
if its incoherent part evolves, as shown in dynamical- 
mean-field theory (DMFT]P 7 for a single-orbital Hub¬ 
bard model on a square lattice in a PM metallic state 
near a Mott transition. Second, Eq. (26) with the ap¬ 


proximate form of g 2 ;acdb{k', 0 ) shows that the intraband 
excitations become dominant compared with the inter¬ 
band excitations. This is because the intraband compo¬ 


nents of Eq. (18) (i.e., a = /?) give larger finite con¬ 
tributions to <7 XX than the interband components (i.e., 
a ^ /?) due to the factor — £* (k) + Fg(fc) in the denomi¬ 
nator of Eq. (18) for w = 0. Third, combining Eqs. (26) 


and (181 with the above second general property, we find 
that cr xx is inversely proportional to the QP damping. 
Note that the dependence of a xx on the QP damping 
can be determined by the dependence of g 2 -,acdb{k', 0 ) since 

A i° 2 ;fea( fc ;°) and A *; 2 -cd (fc; 0) are independent of the QP 
dampin;®. Fourth, due to the CVCs in A^ 2 . ba (k; 0) and 
A x;2 -cd(k; 0), a xx is affected both by the CVC arising from 
the self-energy and by the CVCs arising from the self¬ 
energy and irreducible four-point vertex function, and 
the dominant effect arise from the magnitude changes of 
the currents. This property can be deduced from the 
following arguments: Since A^,° 2 . &a (fc; 0 ) includes the E 


CVC [see Eq. (27)], its effect is the renormalization of 
the group velocity, resulting in a magnitude change of the 
current®!. On the other hand, the effects of the CVCs 
in A x - t 2 - c d(k', 0 ) are not only a magnitude change of the 
current but also an angle change since the CVC arising 
from the irreducible four-point vertex function connects 
the currents at k and fc', which are not always parallel 
or antiparalleP 7 . Those effects on <j xx can be described 


by 

A i; 2 ; 6 a(M)A*; 2 ;c«i(fc; 0 ) 

= l A 2;L( fc )l C0S ^ia( fc )l A 2;cd(^)l COSy> cd (fc) 

H A 2°6a( fc )l C0S < / 3 6a ) ( fc )l A 2;cd( fc )l cos ]> 

(29) 

where lA^^fc)] and |A 2 ;a fc(fc)| represent the magnitudes 

of A i; 2 ;ba( fc ;°) and A *;2 ■cd{k; 0 ), respectively, and (k ) 

and <pab(k) = (k) + A ip a b(k) represent the angles. 

Thus, even for the CVCs in A 2 .. 2 ;c d(fc; 0), the magnitude 
change is dominant for cr xx . 

From those properties, we can deduce the properties 
of the resistivity about the dominant excitations, the de¬ 
pendence on the QP lifetime, and the main effects of the 
CVCs. Since the resistivity is the inverse of the longi¬ 
tudinal conductivity, the dominant excitations and the 
main effects of the CVCs are the same for a XXl and the 
resistivity is inversely proportional to the QP lifetime (in 
the same way for the relaxation-time approximatiorP). 


2. Hall coefficient 


For discussions of the usual Hall effect of a correlated 
multiorbital system for a weak external magnetic field, 
we consider a uniform static external magnetic field along 
the ^-direction, which is so weak that the cyclotron fre¬ 
quency, w c , satisfies w c Tt ran s 1 , and derive an approx¬ 
imate expression of the Hall coefficient in the weak-field 
limit on the basis of the linear-re spons e theory in the 
most-divergent-term approximation 47 -!®!. In this deriva¬ 
tion, we assume that the system has the mirror symme¬ 
tries about the xz- and the yz -plane and the equivalence 
between the x- and the y-directiorPSJ these are valid for 
some 214-type ruthenates without the RuOg distortions. 
Because of the mirro r symmetries and the Onsager re¬ 
ciprocal theorenf 52 ! 53 !, we can treat the Hall coefficient, 
which is generally a third-rank axial tensoi^, as a scalar. 
In addition, because of the equivalence between the x- 
and the y-direction, the Hall coefficient in the linear- 
response theorjE -0 j n the weak-field limit becomes 


Rk = bin ^ 

^xx&yy H—^0 ±1 

= ~Y- bm -jj-. 

(7 2 >0 H 


(30) 


Since we had derived a xx in Sect. II B 1, we need to 
calculate Jim in the linear-response theory with the 

most-divergent-term approximation 4 Al£] in this section. 

To calculate lim rT H ‘, we need to derive the U-linear 
#-*■ o H 

terms of a xy . For that purpose, we use the vector poten¬ 
tial, A, instead of H itself, and derive the q -linear and 

A-linear terms. Thus, the Kubo formula for lim -4r- 

h^o H 
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become^ 47 * 48 ! 


cr. 


1 


lirri ——9 lirn 1 i rn - 

H ^ 0,^0 i(q x A y (q) - q y A x (q))e^ 
s(R)/ 


$xy j (g,w) - $^(9,0) 


(31) 


where (?,(*>) is obtained by <&ffl(q,uj) = 

<h xy (g, ifi n —>■ w + *0+) with 


$ xy (q,iH n )=£ fdr j T dr'e m ^-^ 

x (TrJf ( 9 ,r)jf ( 0 ,r'))// 

= y^ J K xyv (q,iSl n )A v (q). (32) 

V 


and K^l(q,u) is obtained by K^l{q,w) = 
K xyv {q,i^n -> oj + i0+), where K xyv (q,i£l n ) is 
given by 




K xyv (q,ifl n ) = Svajj J dr j drV°” (r ' r) 


(T T J x ( 9 ,r)^^ 


+ J / dr I dr 1 I dT"e l - ln{T - T ‘) 

N Jo Jo Jo 

x (T T J x (q, t)J v ( 0 , T')J v {—q, t")). (34) 


k a,b 




Here J H (0) within the linear response becomes 
J H (0) - J (0) = - ^2^2 A(q) • V k (v k ) ba cl b c ka , (33) 

fc a,fe 


Furthermore, using the three-point vector vertex func¬ 
tion and introducing the irreducible six-point vertex 
functioiPSa we can rewrite Eq. (f34l) as 


j2 T ^ ^ 

K xyv (q,i£l n ) = d,',yjy'52'^2 E Qf : 

k m a.b.A.B U 


d(v k y) 


ba 


G aA (k - , if-m )^x;AB (&- , *6 mi ^+i ^m+n)GBb{k -\-1 ^ra+n) 


+vEEEE G fb (^— i ^m)^-x\ba (^— i ^m j 5 ^m+n)^ad(^+; ^ra+n) 

fe m {a} f,g 

X (^+5 ^ mi > ^m)^cp(^+j ^m)-^-y;dc(^+? ^m+n? ^+5 ^m) 

j2 _^ ^ ^ ^ 

+ 7 vEEE]L G ag{k-\-1 it m+n )A-X]ba (&— j 26 m 5 ^+5 ^m-\-n)G cb(lt— 1 it m ) 
k m {a} f,g 

X (/C-|-, ie m + n , k— , 26 m _| -ri)Gfd{k— I 2€ m _|_ n ) Ay^^/C— 1 2€m+n 1 5 26 m ) 

/ 2T \ 3 

+ (n) E E EE E G Bb{k—i ’itm)-h-x\ba{k — , 26 mo 2€ m _|_ n )(jr a A (&+ > 2€ m _|_ n ) 

k,k',k" m,m',m" {a} {A} f,g,F,G 

X Gqq (&_|_, 2 £ m n ) ^-v\gf (^+ 5 2 € m " 1 7 ^m" ^)Gfp (fc_ , it-m" ^}G Bd{k 5 26 m '_|_ n ) A^;^ c (/c ? 2 £ m '_|_ n ,/c , 26 m /) 

X G c c(k i'i'€m , ')^3^ABCDFGfo+'’^m-\-ni k'—i'i'^mi k 5 26 m ', 1 26 m '_|_ n 5 5 26 m" 1 ^+5 26 m "), ( 35 ) 


with fc+ = k zb The terms of Eq. ([35]) can be given by 
represented by the diagrams shown in Fig. [z] Thus, 

the remaining tasks are to derive the q-linear terms of iv ;o 1 - iin A _ „ a 'i+E^E^E^ 
Kxyviq.itln), to carry out its analytic continuation, and xyi/ ' 2 x y I/,a: N 2^i 2^i 

to combine the result with Eqs. (31) and (32). 


fc m {a} 




A-x;ba(ki ^m+n) ^m+n; ^m) 

Ok y 

2f 1 

26 m ) G ac [(ki itm+n) ■ (*^) 


We first derive the q-linear ter ms of E q. (351 in the We can show the four terms of Eq. (36) as the diagrams 
most-divergent-term approximatiorf 47 ®. As I will ex- in Fig. [3] 

plain in Appendix C in detail, the q-lincar terms are Then, we carry out the analytic continuation of Eq. 
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(a) (b) ( C ) (d) 



FIG. 2: Diagrammatic representations of (a) the first, (b) the second, (c) the third, and (d) the fourth term of Eq. ( |35[ ). x in 
panel (a) represents the momentum derivative. 



FIG. 3: Diagrammatic representation of each term of Eq. (361. The minus signs of the diagrams in panels (b) and (c) are not 


explicitly written, x represents the momentum derivative. 


(36). This procedure can be done in the same way for 
<r vv in Sect. II B 1 since the relevant parameter for ana¬ 
lytic continuations about frequency is the frequency de¬ 
pendence and since the frequency dependence of Eq. (36) 


is the same as K vu (ikl n ) expressed in terms of the three- 
point vector vertex function, 

^ ^ ^ ^ ^ ('Cfcz;)baA t /; C) j(fc, ! 0, i£l n ) 


x G ac (k, ic m-\-n )G d b(k,ie 

m) • 


(37) 


Thus, we obtain A Ki y l(q,u) = K^ y l(q,u_\~ Ki y l(q,0) 
in the most-divergent-term approximation 47 ® within the 
linear order of uj/T: 

AKW(q,u) = - ^(q x 5„, y - q v 5 UiX ) J -^-2 


1 

X — 

N 


¥ 


[ Cj 

A x - t 2-ba{k', 0)^-Aj /; 2;dc(fc; 0) 


k {o} 


G 


(R)/ 




ad ^ dk G ‘ 


it\k) 


(38) 


As described in Sect. II B 1, in the most-divergent- 
term approximation^, the contribution from a retarded- 
retarded or an advanced-advanced pair of two single¬ 
particle Green’s functions is negligible compared with the 
contribution from a retarded-advanced pair. 


Combining Eq. (38) with Eqs. (31) and (132)), we fi¬ 


nally obtain an approximate expression of the dc trans¬ 
verse conductivity in the weak-field limit within the most- 
divergent-term approximation: 


lim 

H-¥ o 


xy 

If 


k ' 


r de ( d/ (e h 

J- oo 2n ^ <9e / 


r ^ 

^ ' |^Aa;j2;fca(^i 0) Ay;2;dc{k] 0) 

{a} V 


x Im 


G™(k)£-G%’(k) 


¥ 


(A), 


dk. 


(39) 


in Sect. II B 1 


Adop ting the similar arguments for a x 
to Eq. (39), we see four general properties for lim <7 jf L . 

First, the QPs near the Fermi level are dominant due to 

the factor (— ). This is the same for a xx . Second, 

the dominance of the intraband excitations also holds 

because of the similar reason for a xx . Note that we can 

obtain the finite intraband components in lim since 

h^o H 


the quantities in the former square bracket in Eq. (39) 
are odd about k x and even about k y due to the com¬ 
bination of the k x derivative in A x - 2 ; ba(k; 0) , and 

the k y derivative in n. y 
ties in the latter are odd about k. 


y derivative in A y -2-dc{k',0), and since the quanti- 


and even about k v 


V 


due to the combination of and a product of the re- 
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tarded and the advanced single-particle Green’s function. 
Third, in contrast to a xx , lim is inversely propor¬ 
tional to the square of the QP damping. This is be¬ 
cause the momentum derivative in a retarded-advanced 
pair leads t o an additional factor of the inverse of the 
QP damping 47 * 48 !. Fourth, the CVCs in A x - 2 -ba{k] 0) and 

A y - 2 -dc(k; 0) affect lim hV, and the dominant effects are 
o H 

an angle change, which is different from the fourth prop¬ 
erty for cr xx . This property arises from the dependence 
of the following quantity on the magnitude and angle 
changes of the currents: 


V 

A X \2\ba(ki 0 ) gj^ Ay2^dc{k^ 0 ) 


= |A 2 ;6a(fc)| COS (fiba(k)\A2-dc(k)\ COS Vdc{k)—^^ 

OKy 

+ |A2;ba(fc)l sin<i ba (fc) d^°| fc ) |A 2; dc(fc)| Sin ip dc (k) 

H A 2 -ba ( k )I cos ^{k)\A 2 . dc (k )I cos ip ( d °J ( k ) 

+ |A 2 ;ba (k )I sin <p { ° a ] (fc) \A 2 .dc{k) | sin ip ( d °J ( k ). 

(40) 


Thus, due to the appearance of dv> g^ k ' > or , the 

angle change of the current causes a more drastic effect 
on lim//_>.o than <r xx . Actually, the importance of 
such drastic effect has been obtained in a single-orbital 
Hubbard model on a square lattic^M 

Combining those properties with the four properties 
for a xx , we can deduce the properties of about the 
dominant excitations, the dependence on the QP lifetime, 
and the main effects of the CVCs. First, the dominant 
excitations are the intraband excitations near the Fermi 
level. Second, the dependence of the numerator and de¬ 
nominator of f?H on the QP lifetime cancels each other 
out in the absence of the band dependence of the QP life¬ 
time, while the cancellation is not perfect in the presence 
of the band dependence. This is because lim#-^ or 
a xx consists of the sum of the corresponding intraband 
components, each of which has the dependence of the 
QP lifetime for the band. Note that the non-perfect can¬ 
cellation is the origin of the temperature dependence of 
i?H of a multiorbital system in the Fermi liquid. Third, 
the main effects of the CVCs on i?n are the magnitude 
change of the current due to A^ . ba {k; 0) in the denomi¬ 
nator of iiu and the angle change of the current due to 
A x - 2 -ba{k\ 0) or Ay- 2 -dc{k-, 0) in the numerator since there 
is the nearly perfect cancellation between the magnitude 
changes due to A x -^ba{k\ 0) and A y ^ 2 -dc{k\ 0) in the nu¬ 
merator and due to the square of A x;2 - c d(k', 0) in the de¬ 
nominator. 


C. FLEX approximation with the E CVC, the MT 
CVC, and the AL CVC 

In this section, after explaining several advantages of 
the FLEX approximation with the CVCs arising from 
the self-energy and irreducible four-point vertex func¬ 
tion, I formulate the FLEX approximation in Matsubara- 
frequency representation for a multiorbital Hubbard 
model in a PM state and derive the CVCs arising from 
the irreducible four-point vertex function in the FLEX 
approximation. In the latter derivation, we first derive 
the irreducible four-point vertex function in Matsubara- 
frequency representation; second, we convert it into a 
real-frequency representation by using the analytic con¬ 
tinuation; third, we calculate part of the kernel of the 
CVCs arising from the irreducible four-point vertex func¬ 
tion; fourth, we derive the Bethe-Salpeter equation for 
the current including the CVCs. Furthermore, I intro¬ 
duce a simplified Bethe-Salpeter equation by approxi¬ 
mating the AL CVC to its main terms. 

To describe the electronic properties near or away from 
a magnetic QCP, I use the FLEX approximation with the 
CVCs arising from the self-energy and irreducible four- 
point vertex function since its following three properties 
are the advantages in describing the electronic trans¬ 
ports. One is that this approximation is one of the con¬ 
serving approxi matio n^SHHI that automatically satisfies 
conservation law d 55 * 56 ) This is powerful to describe trans¬ 
ports since the treatment in keeping conservation laws is 
essential in transports®) Another advantage is that this 
approximation can take account of the many-body ef¬ 
fects due to the self-energy itself and the CVCs arising 
from the se lf-en ergy and the irreducible four-point ver¬ 
tex function 3 10 . In particular, this approximation can 
sufficiently treat the effects of spatial (i.e., momentum- 
dependent) correlation even near a magnetic QCI® toEH. 
Due to this advantage, the FLEX approximation with the 
CVCs can analyze how those many-body effects influence 
the electronic properties beyond random-phase approxi¬ 
mation (RPA), a mean-field-type approximation, and the 
relaxation-time approximation®) where all the CVCs are 
neglectecP, and improve several unrealistic results in the 
RPA; examples of the improvements are a reasonable 
value of U for a magnetic transition and the Curie-Weiss- 
type temperature dependence of the spin susceptibility 
near an AF (QCJI(As described in Sect. II B, the 
CVCs are vital to satisfy conservation law d 3 * 10 * ^ .) The 
other advantage is that the FLEX approximation can suf¬ 
ficiently describe the coherent parts of the single-particle 
Gree n’s functi on for a moderately strong electron corre¬ 
lation 3 10 29 32 ( Actually, the FLEX approximation for a 
single-orbital Hubbard model on a square lattice at U be¬ 
ing a half of the bandwidth is in satisfactory agreement 
with the quantum Monte Carlo calculation about the 
imaginary-time dependence of the single-particle Green’s 
function for several momenta 32 . Although it has been 
proposed in a diagrammatic Monte Carlo calculation— : 
for the same model that diagrammatic expansions based 
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on the Luttinger-Ward functional break down at a large 
U, I believe the above satisfactory agreement 32 remains 
valid since it has been showrP^ that this proposal results 
from an artifact of the technical pathological treatment 
of the noninteracting single-particle Green’s function in 
the diagrammatic Monte Carlo calculation. This suffi¬ 
cient description of the coherent part is very useful to 
analyze the electronic dc transports since, as described 
in Sect. II B, the coherent parts almost dominate the 
electronic dc transports. 

We start to formulate the FLEX approximation for a 
multiorbital Hubbard model in a PM state in a similar 
way for Refs. [33] and [34] A set of the equations in this 
approximation can be obtained by choosing the form of 
the Luttinger-Ward functional as the bubble and the lad¬ 
der diagrams of the multiple electron-hole scattering and 
deriving the effective interaction and the Dyson equa¬ 
tion. First, we can derive the effective interaction in the 
FLEX approximation by considering the bubble-type and 
the ladder-type multiple electron-hole scattering. Since 
we focus on a PM state, it is sufficient to consider the 
following three components: 

Vai= 7j(Ua bcd + u^ bcd ) 

+ 2 E Ua bAB XA BCD {Qii^n)Uc Dcd 

{A} 

~ O E UgbABX < ABCp{ < h ^nWcDcdi (41) 

{A} 

VlU^n) = \{~Ul bcd + U^ bcd ) 

E ^abAB^ABCD^T^nWcDcd 

~ {A} 

^9 E UabABXABCD(Q^^n)U§ Dcd , (42) 

~ {A} 


and 


Xabcd)*?; *Hn) — Xabcd^Qi i^n) 

~ E XabAB{.(l, 'i'^n)U A BC dXcD cd(*}i *^n)- 
{A} 

(46) 


In deriving the effective interaction, we do not need 
to explicitly consider the ladder-type contributions in 
equal-spin-scattering case since those are included in part 
of Eq. (42) as a result of the relation between the 


non-antisymmetrized and the antisymmetrized bare four- 
point vertex function®!. Combining the three compo¬ 
nents and using cr° lS2 , <7° 4S3 , cr SlS2 , and <x S4S3 , we can 
express the effective interaction in the FLEX approxima¬ 
tion as the following single equation: 


ySlS 2 S 3 S4/ Q \ 

V abed (Lh^Ln) 

Uabcd ~ E UabABXABCDiQi^nWcDcd 
{A} 

Uabed +E UabABXABCDiQi ^nWcDcd 

{A} 


2 L 


2 L 


o o 

C7 SiS 2 C7 S 4 S 3 


'SlS2 ^5453 ‘ 


(47) 


Then, using Eq. (471 and excluding the double counting 
of the topologically equivalent term in the self-energy, we 
can derive the Dyson equation, 


G ab (/s,i 6 m ) G ab (k 1 ie m ) 

+ 5 >%(*m £ m )^AB (^b Bb {k 1 Z6 m ), 

A,B 

(48) 


and 

V abcd^ iQ n) = - U*bed 

-EE UabABXABCD(<h ^ n)Uc Dcd > 

{A} s" 

(43) 


where G° fc (fe,«e m ) is the noninteracting single-particle 
Green’s function, 


G° ab (fc,ie m ) = 5^(14°) 


1 


- ea(fc) 


(K^ab, 


(49) 


with 

Xabcdi^Qi ^ ^ G ac {k H - ^m-\-n)G db (^5 ^m)j 

k.m 

(44) 

Xabcd^Q'^n) = XabcdiQ > *Hra) 

+ 5Z XabAB^Q, iQnWABCDXcDediQi 
{At} 

(45) 


and X ab (fc,*e m ) is the self-energy in the FLEX approxi¬ 
mation, 

^ac (k, ie m ) — n EE Vabcd{Qi iQn)Gbd(k 

q,n b,d 

(50) 

with (Ufr)aa, being the unitary matrix to diagonalize 
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e a b{k), and I4 6cd (g, if2„), being 




= -K 


■tttt/ 


abed (?) *^n) ^abed (<Z> *^n) 

— UabABXABCD{q, i^n)Uc Dcd 
{A} 

3 


2 L 

1 r 

"2 . 


Uabed 


UabABXABCD(q ) ^nWc 
{A} 


Ded 


~Uabed + UabABXABCDiQi bXl n )U§ 


Ded 


Ml 


- ^2 U lAbBXABCD{q , 

Ml 


n 

CcDd' 


(51) 


The reasons why the double counting term is the last 
term of Eq. (51) are that the second-order terms in 
V 


tttt 

abed 


(q,iQ n ) and V^ b J d (q,iCl n ) lead to the topologi¬ 
cally equivalent contributions to the self-energy, and that 
Vabcd(q> *^ ra ) contains a relative \ factor arising from the 


coefficient cr+ a 


s l s 2 


Si S o 


sistent set of Eqs. (|44j), (|45 ), (|46j), (|48j), (|50j) , and ( |51[ ) 
with Eq. (49) and the equation to determine /r, 




Solving a selfcon- 


O _ _ OT 1 _ _ ° r 

n * =jf Y.J2 /( £ «( fc )) + [ G aa(k, ie m ) 

k ct km a —1 

- G° aa (k,ie m )], (52) 

we can determine the single-particle or the two-particle 


quantities in the FLEX approximation. Its technical de¬ 
tails for the numerical calculations will be described in 
Appendix D. 

We turn to the Bethe-Salpeter equation for the cur¬ 
rent with the CVCs in the FLEX approximation. The 
derivation consists of four steps. The four steps are to 
derive the irreducible four-point vertex function in the 
FLEX approximation in Matsubara-frequency represen¬ 
tation, to convert it into in real-frequency representation 
by the analytic continuations, to calculate part of the 
kernel of the CVCs, and to combine the part and Eq. 


(28). 


First, we derive the irreducible four-point vertex func¬ 
tion in the FLEX approximation in Matsubara-frequency 
representation. Since the irreducible four-point vertex 
function is generally determined bjtm 

T (1) (k ie k! ie r a ifl 1 - ^ab{k.ie m ) 

1 a6cdv K ’ K > le ™ i 9 ’ lii n) — £Q ■ e r y v 03 ! 

we adopt this equation to the self-energy in the FLEX 
approximation. For the actual calculation, we calculate 
the right hand side of Eq. (53) at q = 0 and fl n = 0, and 


then we label q and fi n so as to represent the electron- 
hole scattering process among an electron of orbital b 
with ( k,ie m ), a hole of orbital d with (k',ie m >), an elec¬ 
tron of orbital a with (k + q,ie m + zfl„), and a hole of 
orbital c with (fc' + q, ie m ' +i£l n ). After the actual calcu¬ 
lation explained in Appendix E, we obtain the irreducible 
four-point vertex function in Matsubara-frequency repre¬ 
sentation in the FLEX approximation: 


with 


and 


^af>ce/(^b k , iCrn 1 ■ q ; i^n) — ^abed k , iCrn 1 : *7: A7 T) ) -)- ^ abcd (fc, ?€ m , k , iCrn 1 j d■ 

+ V ab<!d L2 ( k i ie m , 5 ^n), 


r abed i^m'i — $q, 0 $n,§Vaebd{k k ,26777, 26777,/), 

T 

b i,bcd (k, ie m , k , i€ m '\ q^Sln) = — — y^\ y^] W aBcA . dCbD (q — q , iXl n - n i; —q , — iXl n ') 

q' n> {A} 

X GcA(k + Q , 26777,/+rc/ ^)GrBD {k H - Q ■> ^ra+rc/)) 


r i 1 b ) < y L2 (k,ie m ,k',ie m r } q,in n ) = - ^ 5Z W aBAd-CcbD(-q^ q - q\ ~i^n+n') 

q' n' {A} 

X GAc{k Q 1 iCm' —n' ^)GbD {k Q Q 1 ^m+n+n')? 


(54) 

(55) 

(56) 


(57) 


r 


where W^ d . ABCD {q 1 ,in ni -,q 2 ,iCt n2 ) is 


W abed; ABCD (<? 1 > 5 <? 2 , i^n 2 ) 

= \ ftl bcd {qi , itt ni ) ft ABC D (92, itt n2 ) 

+ ^ftabed(q^ in n 1 )ftABCD(q2,i^n 2 ) 

— u n u n 

u acbd u ACBD1 


with 

ftabcdiq'i i^n') = U^ bcd + ^ U^ bCD XcDAB(q'i ^n')U A BcdJ 

{A} 

(59) 


(58) 
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(a) MT term 


(b) AL1 term 


(c) AL2 term 


and 



fc'; 0) = F^ 2 (e, e'; T) (^)^EE 

q' {A} 

/ oo 

dw , i^ L2 (e, C ',a/;T) 

-oo 


' —oo 
^AL(RA) 


FIG. 4: Diagrammatic representations of the irreducible four- 
point vertex functions in the FLEX approximation. Neglect¬ 
ing orbital indices and relabeling momentum and frequency 
variables, we can show that each diagram is equivalent to the 
corresponding diagram in Ref. 1371 


and 

NabcdWi i^ri) = Uab cd ~ E ^CdXgDAB) G AB cd- 

{A} 

(60) 


x&h W) 
x ImG Ac( fc/ - 9') ImG E2( fc + v')> 

(64) 

where F™ T (e,e'-,T), F^\e,e'-,T), F t f L 1 (e, e', w'; T), 
F c ^ L 2 (e, e'; T), and F t ^ L 2 (e, e', u/; T) are, respectively, 


F c I f T (e,e , ;T)=(coth 

F c E(e,e , ;T)=fcoth 


e — e 
2 T 
e'-e 


tanh ^)’ 


ct 
dALI 


We can represent the terms of Eqs. (55), (56), and (57) 
as the diagrams of Figs. [4|a),[4](b), andQc), respectively. 


Second, we carry out the analytic continuations of 
Eqs. (55), (56), and (57) to convert these into real- 


frequency representation. This is because the irreducible 
four-point vertex functions in real-frequency representa¬ 
tion are necessary to calculate part of the kernel of the 

CVCs, J-zll-dcni^- k'\ 0) [see Eq. (J2sj)] - Carrying out the 
analytic continuations, we obtain the MT, the AL1, and 
the AL2 terms for regions 22-11, 22-III, and 22-IV (see 
Appendix F). 

Third, using the MT, the AL1, and the AL2 terms 
in regions 22-11, 22-III, and 22-IV, we can calculate 
Elide k'\ 0) in the FLEX approximation. Since the 
irreducible four-point vertex function is the sum of the 
MT, the AL1, and the AL2 term, ElidCD^’ k'\ 0) i n the 
FLEX approximation is given by 

E!;Ld(*b k'; 0 ) =j£>%£(k, k'- 0 ) + E 1 E 1 (k, k 0 ) 

+ J$a^(k,k'-, 0 ), (61) 

with 

ElSc T d(fc. k'-, 0) = F c r (e, e'; T)2ilmV™ d {k - k'), 

(62) 
N 

q' {A} 

/ OO 

du'F£ L1 (e,e,u';T) 

-OO 

X ImG^l(k' + </)ImG^(fc + q 1 ), 

(63) 


k'; 0) = F^(e, e'; T) (^)^EE 


(e,«W;T) = tanh 


F 


-AL2 


(e, e'; T) =^coth 6 — tanh 


— tanh — 

2 T 2 T 

j' + e , u/ + e' 

- tanh ^“ 

2T ) 1 


2 T 


2 T 


and 


-iAL2 


(e, e 7 , u/; T) = (tanh ■ 


2 F 


e w — e 

— tanh 


2 F 


(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


In Eqs. (62) (64), we have used the relations of the effec¬ 


tive interaction and the single-particle Green’s functions 
due to the time-reversal and the even-parity symmetry; 
in more general case, we should not use the relations such 
as V^ d (k - k') - V^ d (k - k') = 2HmV™ d (k k') and 

G B£>( fc + ?') - G E( fc + q’) = 2 ilm G { B l(k + q '), and 
should retain the differences between the retarded and 
the advanced quantities. 


Fourth, substituting Eqs. (61) with Eqs. (62), (63), 


and (64) into Eq. (28), we obtain the following Bethe- 


Salpeter equation with the CVCs in the FLEX approxi¬ 
mation: 

A , ; 2 ;cd(fc; 0 ) =A g. cd (fc; 0 ) + AA^ cd (fc; 0 ) 

+ AA^ 2 lc d {k; 0) + AA^ ; 2 f c d(fc; 0), (70) 

where AA ™J. cd (k;0) is the MT CVC, 

AA“E( fc ;°) =wEE / b^- F c 1 f r (ce / ;2 1 ) 

X Im EcE( fc - k')KncD{k'-, 0 ), (71) 


AA^. 2 1 cd (fc; 0) is part of the AL CVC, 


'; 2 \cd 
i AL1 


i 




noo /*oo 

de' / duj' 


xff( e / ; r)C( £ ,^';r) 

x ImG^E {k' + </)ImG^E,(^ + 
xi ; 2 ;CD (fc';0), (72) 
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and AA^. 2 2 cd (fc; 0) is the other part of the AL CYC, with 


aa ^; 2L( fc ; 0) = - 


1 


4t r 2 V 2 


EEE 


de' 


duj' 


1 — OO 
J . 


k' ,q' {A} {A'> 

xFf( £ , e ';T)C(e,^';T) 

xW^A'D]ccdD'(-Q'^-Q') 
x ImG^p, ( k 1 — (k + q ') 

xY ; 2 ;CD (fc';0), (73) 


with 


A„ ; 2 ;Cd(A;'; o) = ^ 92-CABD{k' ] 0)A l y-2-AB(k'- 1 0). (74) 


A,B 


Equations (fTlj) ([73]) show that the MT and the AL CVC 
connect the currents at different momenta; for example, 
the MT CVC connects the current at k with the current 
at k'. 

In the actual numerical calculations, instead of the 
above Bethe-Salpeter equation, I use the simplified 
Bethe-Salpeter equation where the AL CVC is simplified 
by only its main terms. The main terms of the AL CVC 
can be determined by using the following two properties 
satisfied in the present model: The terms arising from 
U are dominant compared with the terms arising from 
the other interactions in a realistic parameter set (i.e., 
U > U ', U > Jh, and U > J'); the intraorbital compo¬ 
nents of the current are larger than the interorbital ones 
due to the large intraorbital hopping integrals compared 
with the interorbital hopping integrals (i.e., the larger 
intraorbital components of the group velocity). Namely, 
the main terms of the AL CVC are given by the sum of 
the following two quantities: 


AA 


ALl 

is;2’,cd 


(fc;0) 


and 


~6 C 


1 


’ 47t 2 V 2 


E 


/*oo /*oo 


de' / dJ 


' —OO J —OO 


k' ,q' * 

x IV c AL(RA) (- 9 ';-g') 
x ImG^(fc / + g')ImG^(fc + q’) 
xi ; 2 ;cc (fc';0), (75) 




1 

l 4^N 2 


E 


rOO /»oo 

/ de' / dJ 

OO J— OO 


k' ,q' * 

xF AL 2 (e, e ';T) J F AL 2 (e, e ', W ';T) 
x W^ HKA \-q'-,-q') 
x ImG^(fc' - q)lmG^\k + q') 
x A„ ; 2 ;cc (fc'; 0), (76) 


where IV C AL(RA ^(—g'; —g') is given by 

W AHKA) { _ q/ ._ q/) = 3 -N?£(- q >)N? c £\- q ') 

+ (-q)Ngif i-q') - U 2 , (77) 


E S cS } (-90 = U + U\ s ™ (-g0, (78) 

(-90 = U- U 2 x^ (-gO- (79) 


More precisely, by using the former of the above 
two properties (corresponding to considering 
only the terms arising from U), we can replace 

WcB^CA^DC'dD'i-d 1 '’ -90 of the ALl term and 
W^B^A’D-C'CdD ’(~90 — 90 AL 2 term by, respec¬ 
tively, S b ^ c S c ,cSa' , c SD,dSc ,dSD' ■ -q') 

and Sb\ c 5a\c&d,c8c' 4&c,ddD' ,dW^ctuld (—gO —90; 

furthermore, using the latter property, we obtain Eqs. 
[75]) and [76]). Solving Eqs. [70]), [7l]), [75]), and [76]) 
with Eq. (271 self-consistently, we can determine the 
current including the CVCs arising from the self-energy 
and irreducible four-point vertex function in the FLEX 
approximation. I will describe the technical remarks to 
numerically solve those equations in Appendix G. 


III. RESULTS 

In this section, I show the results of the magnetic prop¬ 
erties, the electronic structure, and the transport prop¬ 
erties for a PM state of the multiorbital Hubbard model 
away from or near the AF QCP. In Sec. Ill A, I present 
the results of the magnetic properties in the FLEX ap¬ 
proximation. From those results, we discuss the dom¬ 
inant fluctuations, the static and the dynamic proper¬ 
ties of the spin susceptibility, the role of each t 2 g orbital, 
and the effects of the spin fluctuations on the imaginary 
part of the retarded effective interaction. In Sec. Ill B, 
to discuss the effects of the self-energy on the electronic 
structure, I show the results of the FS, the mass enhance¬ 
ment factor, the unrenormalized QP damping, and the 
QP damping in the FLEX approximation. In Sec. Ill C, 
we discuss the main effects of the AL CVC on the inplane 
resistivity, p a b , and the Hall coefficient in the weak-fielcl 
limit, i?H, in the FLEX approximation with the X CVC, 
the MT CVC, and the main terms of the AL CVC and 
more simplified three cases. In addition to the temper¬ 
ature dependence of those transport coefficients, I show 
the orbital depedences of a xx and lim j n order to 

determine the role of each t 2 g orbital. 

I obtained the results of this section by the numer¬ 
ical calculations using the techniques explained in Ap¬ 
pendixes D and G and converting the quantities obtained 
in the FLEX approximation in Matsubara-frequency 
representation to the corresponding quantities in real- 
frequency representation by the Pade approximation 61 62 . 
In the numerical calculations, I used N = 64 x 64 meshes 
and M = 1024 Matsubara frequencies, set Aej = Ac), = 
A J-„ = 0.0025 eV, e c = 0.2 eV, J' = J H , U’ = U - 2 J H , 
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and chose Jr, U, and T as parameters. (A e g , Ac',, and 
Auij,, are the intervals of the discretized real-frequency 
integrals, and e c is the cut-off frequency in the discretized 
real-frequency integrals.) The parameters of Jr, U, and 
T were chosen as follows: I put ,/r = F except the anal¬ 
ysis of the dominant fluctuation, in which the value of Jr 
was chosen in the range of 0 < Jr < 5; I considered the 
case of U = 1.8 or 2.1 eV as, respectively, cas e! 10 * 38 ! away 
from or near the AF QCP except the results in noninter¬ 
acting case; I considered several values of T in the range 
of 0.006eV < T < 0.03eV. In addition, in the conver¬ 
sion by the Pade ap prox imation, I numerically solved its 
recursive procedur e ! 61 * 62 l using the quantities at the low¬ 
est four Matsubara frequencies; for example, we obtained 
E^'(fc,e) by adopting that recursive procedure to a set 
of Y, ab (k , ie m ) at m = 0,1, 2, 3 in the FLEX approxima¬ 
tion. Note that the advanced quantities are obtained by 
using the relations such as £^(fc,e) = £ a R ^(fc, e)* due 
to the time-reversal and the even-parity symmetries. 


A. Magnetic properties 

In this section, I show four main results about the 
magnetic properties. First, the dominant fluctuations 
are the spin fluctuations. Second, an increase of elec¬ 
tron correlation leads to the enhancement of low-energy 
spin fluctuation at q = Qic-AF = (§j7r, M 71 ")- Third, the 

diagonal and the nondiagonal component of X^abb (?) 
q = Qic-af contribute to the enhancement of the spin 
fluctuation at q = Q ic-af, and the diagonal component 
of the d xy orbital is largest. Fourth, the orbital depen¬ 
dence of the effective interaction is determined by the 
orbital dependence of the spin fluctuation. 

We first determine the dominant fluctuations in the 
present model. For that purpose, we analyze the effects of 
electron correlation on A^ ax (q ) -1 and A^ ax (q) _1 , the in¬ 
verses of the maximum eigenvalues ! 63 1 64 1 of X a 6 cd (g, 0 ) and 
Xaf>cd(9> 0), respectively. This is because by analyzing the 
dependence of A^ lax (qr ) —1 and A^ ax (qr ) —1 on U and J H , 
we can determine the dominant fluctuations among four 
kinds of fluctuations, i.e. charge fluctuations, spin fluc¬ 
tuations, orbital fluctuations, and spin-orbital-combined 
fluctuation!^! (for more details see Appendix H). I 
show A^ ax (q ) _1 and A^ ax (< 7) _1 at T = 0.006 eV and 
U = 1.8 eV for several values of Jr in Figs. [5](a) and[5](b), 
respectively. We see that as Jr increases, A^Jq )' 1 
monotonically decreases, and A^ ax (< 7) _1 monotonically 
increases. This behavior is characteristic of the enhance¬ 
ment of spin fluc tuatio ns and the suppression of the 
charge fluctuation ! 65 * 66 ! (see Appendix H). The similar 
results are obtained at U = 2.1 eV, as shown in Figs. |5](c) 
an did). Since approaching the inverse of the maximum 
eigenvalue towards zero characterizes the enhancement of 
the susceptibility, the results in Figs. § a )-§d) show that 
spin fluctuations are dominant at U = 1.8 and 2.1 eV in 
the present model. This can be understood by consider¬ 


ing the following three facts: the noninteracting suscep¬ 
tibility for a = c and b = d, xi&a&(*Z> *^n), becomes very 
large in the present model since G'S (k. icu m ) is larger 
than G^(k,icj m ) for b yf a due to the large intraorbital 
hopping integrals compared with the interorbital ones; 
the interactions between the different kinds of fluctua¬ 
tions may be generally very weak in the FLEX approx¬ 
imation due to lack of the vertex corrections of the sus¬ 
ceptibilities; the terms arising from U cause the strongest 
enhancement of the susceptibilities. Namely, due to those 
facts, the intraorbital components of XabcdiQ i*^n), be. 
X aaaQ (q, zfl n ), are strongly enhanced, resulting in the 
larger enhancement of spin fluctuations than the other 
fluctuations. Hereafter, we fix the value of Jr at Jr = 
Then, for a deeper understanding of spin fluctuations 
in the present model, I analyze the static and the dy¬ 
namic properties of the spin susceptibility as a function 

of w, x S(R) (<Lw) = Ea.bXaiwk'Lw)- For the analysis of 
the static property, I show the momentum dependence of 
X S (<7, 0) = x s ( R )(g, 0) at T = 0.006 eV for U = 1.8 and 
2.1 eV in Fig. | 6 ](a). The result shows that as U increases, 
the spin fluctuation at q = Q ic-af is most strongly 
enhanced and the enhancement at q = ( 0 , 0 ) is much 
weaker. That strongest enhancement can be understood 
as the combination of the merging of the nesting vectors 
of the d xz /y Z and d xy orbitals around q = Qic-AF due to 
the mode-mode coupling for the spin fluctuations around 
q = Qic-AF and the nesting instability at q = Qic-af 
due to the RPA-type scattering process, as explained in 
Ref. (38l Next, for the analysis of the dynamic prop¬ 
erty, I show the frequency dependence of Imx s ( R H<b w ) 
for several values of q at T = 0.006 eV for U = 1.8 and 
2.1 eV in Figs. §b) and[ 6 ](c). These figures show that 
low-energy spin fluctuation at q = Qic-af is dominant 
in the dynamic properties at U = 1.8 and 2.1 eV, and 
that the intensity at q = ( 0 , 0 ) is very small. 


Moreover, I analyze the role of each t 2 g orbital in dis¬ 
cussing the spin fluctuations. Figures[7](a)-[7](d) show the 

frequency dependences of Imy^ (Qic-af, w) at J H = § 
for (T,U) = (0.006,1.8), (0.006,2.1), (0.02,1.8), and 
(0.02,2.1) (eV). We see that not only the diagonal but 
also the non-diagonal components are enhanced, and that 
the largest component is the diagonal one of the d xy or¬ 
bital. First, the enhancement of the diagonal components 
arises from the combination of the large diagonal com¬ 
ponents of the noninteracting susceptibility of the t 2 g or¬ 
bitals around q = Q ic-af, the merging of the nesting vec- 
and the d xy orbital around q = Qic-af> 


tors of the d x 


■-/v z 


and the larger enhancement due to the terms arising from 
U than the other terms. Next, the nondiagonal compo¬ 
nents are enhanced due to the terms including Jr and 
the diagonal components since XaabbiQ’ ^n) for a ^ b are 
enhanced mainly through Xaaaa{q,i^n)Uxl abb {qGU n ) + 


Xa 


i{q, *H n ) Juxfbbbiq^ *A1„) [see the second term of Eq. 


(45)]. Then, the diagonal component of the d xy orbital 
becomes largest due to the following three properties: 
the diagonal components of the noninteracting suscepti- 
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(a) (b) (c) (d) 



FIG. 5: Momentum dependence of A| lax (q) 1 and An lax (g) 1 at T = 0.006 eV and U = 1.8 and 2.1 eV for several values of 
Jh. 



FIG. 6: (a) Momentum dependence of x s (g,0) = x s l R l(g, 0) at T = 0.006 eV, U = 0, 1.8, and 2.1 eV and Jh = ^f, and 

frequency dependence of Imx s(R ' ) (q, uj) for q = Qic-af, Q[ c _af = (n, §§7r), and Qfm = (0, 0) at T = 0.006 eV and J H = ^ for 
(b) U = 1.8 eV and (c) U = 2.1 eV. 
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FIG. 7: Frequency dependence of Imx a ^(Qic-AF, w) 

at J H = f for (T, 17) =(a)(0.006,1.8), (b)(0.006,2.1), 
(c)(0.02,1.8), and (d)(0.02, 2.1) (eV). 


FIG. 8: Frequency dependence of ImV^^^Qic-AF,w) 

at J H • f for (T, U) =(a)(0.006,1.8), (b)(0.006,2.1), 
(c)(0.02,1.8), and (d)(0.02, 2.1) (eV). 


bility are larger than the non-diagonal components due 
to the large intraorbital hopping integrals; the noninter¬ 
acting susceptibility of the d xy orbital is larger than that 
of the d xz / yz orbital due to the larger DOS 10 of the d xy 
orbital; the enhancement due to the terms arising from 
U is largest in the terms arising from the Hubbard inter¬ 
action terms. 

Finally, we see the effect of the spin fluctuations on 
the imaginary part of the retarded effective interaction 


of the FLEX approximation. The reason why that effect 
is analyzed is that its understanding is useful to under¬ 
stand the effect of the spin fluctuations on the MT CVC 
since the imaginary part of the retarded effective interac¬ 
tion is part of the kernel of the MT CVC [see Eq. (p7T|)] . 
For that analysis, it is sufficient to present ImV oaf ; f) (q, ui) 
since the other orbital components are much less im¬ 
portant. This is due to the facts that the dominant 
fluctuations are the spin fluctuations and that their 
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FIG. 9: FSs at T = 0.006 eV and J H = f for (a) U = 1.8 
eV and (b) U = 2.1 eV. In panels (a) and (b), the FS sheets 
at U = 0 eV are shown by the dotted lines for comparison. 


contributions to the effective interaction, V ac bd(Q , i^r, 
are given by E A,c U acAAXAACC^ in n)Ul Cbd _ 


fia,c3b,d Ej4,C ^aaAAX^iACci^’ [see Eq. (51)]. 

Figures [ 8 ja)-[sj]ci) show the frequency dependence of 

Im Ea R bi(QiC-AF,w) for (T, U) = (0.006,1.8), (0.006,2.1), 
(0.02,1.8), and (0.02,2.1) (eV). The obtained orbital de¬ 
pendence is similar to that for Imx®^(q, w). Thus, the 
spin fluctuations lead to the main contributions to the 
MT CVC in the present model, and the orbital depen¬ 
dence of the MT CVC is determined by the orbital de¬ 
pendence of the spin fluctuations. 


B. Electronic structure 

In this section, I show four main results about the elec¬ 
tronic structure. First, the topology of the FS remains 
the same as the noninteracting one even including the 
FS deformation due to the real part of the self-energy 
in the FLEX approximation. Second, the mass enhance¬ 
ment of the d xy orbital is larger than that of the d xz / yz 
orbital in a wide region of the parameter space in the 
present model. Third, the unrenormalized QP damping 
of the d xy orbital becomes larger than that of the d xz / yz 
orbital. Fourth, the orbital dependence of the QP damp¬ 
ing is mainly determined by the orbital dependence of 
the unrenormalized QP damping. 

I begin with the effects of the real part of the self¬ 
energy in the FLEX approximation on the FS and the 
mass enhancement factor. I determine the FS by cliag- 
onalizing [e a b(k) + ReE^j/ (fc, 0 )], where p in e a b(k) has 
been determined by Eq. ( [52| . 

First, we see from Figs. da) and d^) how the FS 
is modified with increasing U. Those figures show that 
the modification is slight. Thus, the real part of the 
self-energy in the FLEX approximation does not change 
the topology of the FS sheets (i.e., whether each sheet is 
electron-like or hole-like). This result can be understood 
by considering two facts that the occupation numbers of 
the d xz / yz and the d xy orbital do not become very close 
to integers, and that the van Hove singularity of the d xy 
orbital does not cross over the Fermi level. Note, first, 
that the occupation numbers of the d xz / yz and the d xy 
orbitals are 1.36 and 1.28, respectively, at U = 1.8 and 


2.1 eV; second, that if the van Hove singularity crosses 
over the Fermi level, the 7 sheet touches the boundary of 
the Brillouin zone at k = (77 0) or (0, 7 r). 

Next, we show the mass enhancement factor, 

Zaik)- 1 = 1 - OReS £ ,(fc ’") |^o, at U = 1.8 and 2.1 eV 
in Figs. From those figures, we find three 

properties about the orbital, temperature, and momen¬ 
tum dependences of z Q (fc) -1 . The first property is that 
the mass enhancement of the d xy orbital is always larger 
than that of the d xz / yz orbital for all the temperatures 
considered. This arises from the stronger spin fluctua¬ 
tions of the d xy orbital than those of the d xz / yz orbital, 
as explained in Ref. 1381 Combining this result with the 
similar orbital dependences of z a (k)~ 1 as a function of 
Jh (in 0 < Jh < y at T = 0.006 eV and U = 1.8 
eV), we deduce that the larger mass enhancement of the 
d xy orbital is realized in a wide region of the parame¬ 
ter space of the present model for a PM state in the 
FLEX approximation. It should be noted that although 
the spin fluctuations of the d xy orbital enhance z a (k) _1 
of not only the d xy orbital but also the d xz / yz orbital, 
the enhancement for the d xy is larger in a realistic set 
of the Hubbard interaction terms. This is because the 
spin fluctuations of an orbital cause the enhancement of 
z a (k)~ 1 of the orbital proportional to the U 2 terms of 

I Ea,B UaaAAXAABB^Q’ ^n)U% Baa in Vaaaa in 

V a a(fc,fom), and the enhancement of z a (k)~ l of another 
orbital proportional to the J B te rms. Then, the second 
property found in Figs. [l0](a) Jl0](d) is that the tempera¬ 
ture dependence is weak other than the case for the d xy 
orbital at U = 2.1 eV. This results from the more sig¬ 
nificant enhancement of the spin fluctuations of the d xy 
orbital with decreasing temperature [see Figs. 0 a)- 0 d)], 
and suggests that the mass enhancement of the d xy or¬ 
bital may remain larger even at lower temperatures than 
the temperatures considered. The third property of Figs. 
[To];a)-[To]: d ) is that the momentum dependence is negli¬ 
gible for the d xz / yz orbital, while the d xy orbital has the 
weak momentum dependence. This is due to the differ¬ 
ence between the quasi-one-dimensionality of the d xz / yz 
orbital and the quasi-two-dimensionality of the d xy or¬ 
bital: only the d xy orbital has the van Hove singularity 
due to the saddle points at k ss (|| 7 T, 0 ) and ( 0 , §§tt), 
resulting in a larger mass enhancement 67 . Since this re¬ 
sult shows that the momentum dependence of the mass 
enhancement factor is not important to discuss the mag¬ 
nitude difference of the mass enhancement, the present 
analysis is sufficient for that discussion. 


Then, we turn to the effects of the imaginary of the 
self-energy on the unrenormalized QP damping, 7 a (k) = 
—Im£aa(fe, 0). From the results shown in Figs. 11(a)- 
|TTj(d) , we see three main features. The first one is about 
the orbital dependence: the magnitude for the d xy orbital 
is about three times as large as that for the d xz / yz orbital. 
This arises mainly from the larger DOS and stronger spin 
fluctuations of the d xy orbital. Note, first, that a ratio 
of the noninteracting DOSs of the d xy and the d xz / yz or- 
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(a) (b) (c) (d) 



eV, (b) a = d xy and U = 1.8 eV, (c) a = d xz and U = 2.1 eV, and (d) a = d xy and U = 2.1 eV. Zd yz (k) 1 is given by 

Z d yz {kx, k y ) =z d xz (ky,k x ) 




FIG. 11: Temperature dependence of 7 a (fc) = — lm£ia'(fc, 0 ) for several k at Jh = ^ for (a) a = d xz and U = 1.8 eV, (b) 
a = d xy and U = 1.8 eV, (c) a = d xz and U = 2.1 eV, and (d) a = d xy and U = 2.1 eV. The inset in panel (a) or (b) shows 
7 a (k) against T 2 below T = 0.01 eV, and the inset in panel (c) or (d) shows "f a {k) against T 0 5 below T = 0.01 eV. 7 d y ,{k) is 
given by 7 ^ (k x , k y ) = j dzz ( k v , k x ). 


bitals on the Fermi level is about 2. second, that due 
to the similar reasons for 2 a (fc) _1 , the spin fluctuations 
of the d xy orbital cause a larger enhancement of 7 Q (fc) 
of the d xy orbital in a realistic set of the Hubbard in¬ 
teraction terms. The second main feature is about the 
temperature dependence: the unrenormalized QP damp¬ 
ings of the d xz / yz orbital at U = 1.8 eV show the T 2 
dependence at low temperatures; the T 0 ' 5 dependence of 
7 a (fc) for the d xz / yz orbital is realized for k = |1 tt) 

at U = 2.1 eV; the unrenormalized QP damping of the 
d xy orbital at k = (| 7 r, 0 ) is proportional to T linear at 
U = 1.8 and 2.1 eV. The T 2 dependence is due to the for¬ 
mation of long-lived QPfP^ the T °' 5 dependence results 
from the hot-spot structured due to the enhanced AF 
spin fluctuation, as explained in Ref. ID; the T-linear 
behavior emerges as a result of the existence of the van 
Hove singularitjd The third main feature is about the 
momentum dependence: the unrenormalized QP damp¬ 
ing of the d xy orbital depends weakly on momentum; the 
momentum dependence for the d xz / yz orbital is negligi¬ 
ble. This arises from the considerable difference in the 
momentum dependence of the single-particle spectrum 
function due to the existence of the van Hove singularity 
only for the d xy orbital. 

Finally, we analyze the effects of the combination of the 
real and the imaginary part of the self-energy on the QP 
damping, 7 *(fc) = z a ik)"la{k). From the results shown 
in Figs. we see that even for the QP damp¬ 

ing, the larger magnitude for the d xy orbital is realized. 
This is due to the larger difference in the unrenormalized 
QP damping compared with the difference in the mass 


enhancement factor, and suggests that the QPs of the 
d xz / yz orbital are more coherent than the QPs of the d xy 
orbital in the present model. In addition, we find the T 2 
dependence for the d xz / yz orbital at low temperatures at 
U = 1.8 eV, the deviation from the T 2 dependence for 
the d xz / yz orbital at U = 2.1 eV, and the similar mo¬ 
mentum dependence of the QP damping to that of the 
unrenormalized QP damping. 


C. Transport properties 


In this section, I show three main results about the 
transport properties. First, the main results in the previ¬ 
ous studies?--^ remain qualitatively unchanged even in¬ 
cluding the main terms of the AL CVC. Second, the tem¬ 
perature dependences of p a b and I?h near the AF QCP 
consist of two regions, high-temperature region, where 
only the E CVC is sufficient, and low-temperature re¬ 
gion, where only both the E CVC and the MT CVC 
are sufficient. Third, in contrast to the case near the 
AF QCP, the effects of the MT CVC on p a b and Ru at 
low temperatures are different in case away from the AF 
QCP: only for I?h, the effects are considerable. 

To analyze the main effects of the AL CVC on p a b 
and .Rh, we consider four cases, named MT+AL CVC 
case, MT CVC case, E CVC case, and No CVC case. 
In the MT+AL CVC case, we take account of the E 
CVC, the MT CVC, and the main terms of the AL 
CVC: A„ ; 2 ;c£ z(fc; 0) in Eq. (701 includes those CVCs, and 
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(a) (b) (c) 



(d) 



FIG. 12: Temperature dependence of 7 *(fc) = z a (k)"f a (k) for several k at Jh = ^ for (a) a = d xz and U = 1.8 eV, (b) a = d xy 
and U = 1.8 eV, (c) a = d xz and U = 2.1 eV, and (d) a = d xy and U = 2.1 eV. 7 a(k) = T are shown by the dotted lines to 
discuss whether the QP damping is cold-spot-type or hot-spot-type, and the insets show 7 *(fe) against T 2 below T = 0.01 eV. 
7 d yz ( k) is given by 7 ^ {kx, fey) = 7 d xx (*+ k *)■ 



T(eV) 


FIG. 13: Temperature dependence of p ab at U = 1.8 eV and 
Jh = jr in the four cases. The inset shows p a b against T 2 
below T = 0.01 eV. 


A^. o 6 (fc; 0) in Eq. (E7J) includes the E CVC. In the MT 
CVC, we neglect only the AL CVC and take account 
of the other CVCs: the change from the MT+AL CVC 
case is neglecting the AL CVC in A„ ; 2 ;c d(fc; 0). In the E 
CVC case, we take account of only the E CVC among the 
CVCs: A l/; 2 ;cd(fe; 0) becomes the same as A^, a 6 (fc; 0). In 
the No CVC case, we neglect all the CVCs: A 1 /; 2 ; C d(fc; 0) 
and A^ a b(k', 0) are determined only by the noninteract¬ 
ing group velocity. 


1. In-plane resistivity 

We begin with p a b{= &xx = a yy) awa Y from the AF 
QCP. We show the temperature dependence of p ab at 
U = 1.8 eV in the four cases in Fig. [T3j and find 
two main features. One is that the T 2 dependence be¬ 
low T = 0.008 eV holds even in the MT+AL CVC 
case. This can be understood that the CVCs little af¬ 
fect the power of the temperature dependence of the re¬ 
sistivity. This is because the main effects of the CVCs 
on the resistivity arise from the magnitude changes of 
the current [see Sect. II B 1] and because the magni¬ 
tude changes appear in the equation of the resistivity 


as 


|A^ d (fc)|+A|A 2;cd (fe)| ~ |A<°> d (fc)| 


( 1 - 


A|A 2;cd (fc)| 


), where 


is not large. The other main feature is that the 

value of p ab in the MT+AL CVC case becomes smaller 
than that in the MT CVC case and nearly the same as 
that in the E CVC case. This is due to the small ef¬ 
fects of the MT and the AL CVC; for high-temperature 
region, the small effects arise from the dominance of the 
QP damping compared with the spin susceptibility in de¬ 
termining the kernels of those CVCs; for low-temperature 
region, the small effects arise from the combination of the 
not large spin susceptibility and the partial cancellation 
between the effects of the MT and the AL CVC. The 
more detailed explanations about those are as follows: 
In discussing the effects of the MT and the AL CVC, 
the relative values of the spin susceptibility and the QP 
damping are important since the kernels of the MT and 
the AL CVC contain the spin susceptibility and the in¬ 
verse of the QP damping [see Eqs. (701 (73)]. Due to 
this property, at high temperatures, the kernels become 
small since the QP damping is large; thus, the effects of 
the MT and the AL CVC are small for high-temperature 
region. Furthermore, although the effects of the MT and 
the AL CVC are separately non-negligible at low temper¬ 
atures since with decreasing temperature the QP damp¬ 
ing decreases and the spin susceptibility remains almost 
unchanged, the effects of the AL CVC reduce the effects 
of the MT CVC as a result of the difference in the momen¬ 
tum dependence; due to this reduction, the total effects 
of the MT and the AL CVC are small. Such property 
due to the difference in the momentum dependence can 
be easily seen from a simple and sufficient case of the 
second-order perturbation theory for a single-orbital sys¬ 
tem since the momentum structure of each diagram of 
the MT, AL1, and AL2 terms remains the same as in the 
FLEX approximation: in this case, the MT CVC is given 
by Yhk' q Ao(fe, k'\k'+q, k-q)& k ^ q (e), and the AL1 and 
AL2 CVCs are J2 k , q A 0 (fc, k'\k' + q,k- q)[<*> kl+q (e) - 
3>fc/(e)] (for more details, see Ref. [15]) ; since $fc(e) is odd 
about momentum, the difference in the sign of q leads to 
the partial cancellation of the effects of the MT and the 
AL CVC. 

We next discuss the role of each f 2s orbital in determin- 
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(a) (b) (c) 



FIG. 14: Temperature dependence of (a) a xx and orbital-decomposed components in the MT+AL CVC case at U = 1.8 eV 
and the orbital-decomposed components of (b) the d xz and d yz orbitals and (c) d xy orbital in the four cases at U = 1.8 eV. 



FIG. 15: Temperature dependence of p a b at U = 2.1 eV and 
Jh = ^ in the four cases. 


ing p ab away from the AF QCP. For that purpose, I show 
the orbital-decomposed components of a xx , the d xz + d yz 
component and the d xy component; the former is ob¬ 
tained by replacing Xq a }=i in Ec l- (26) by Xq a }=i’ and 
the latter is obtained by replacing X]{a}=i hi Ec l- (26) 
by S{o}= 3 - As explained in Ref. TUJ only those com¬ 
ponents are sufficient in the present model since those 
(diagonal) components are larger than the non-diagonal 
ones due to the large intraorbital hopping integrals com¬ 
pared with the interorbital hopping integrals. First, we 
see from Fig. 14(a) that <J xx / yy is determined almost by 
the component of the d xz / yz orbital, and that the contri¬ 
butions from the component of the d xy orbital are very 
small. Namely, the d xz + d yz component remains domi¬ 
nant even with the main terms of the AL CVC. We also 
see from Figs. E3b) and [l4|c) that the values in the 
MT+AL CVC case are nearly the same in the E CVC 
case. Thus, the E CVC is sufficient for discussions about 
the orbital dependence away from the AF QCP. 

From the results at U = 1.8 eV, we deduce, first, that 
the main results obtained in the previous studie^SEII 
away from the AF QCP, the T 2 dependence of p ab at 
low temperatures and the dominance of the d xz / yz or¬ 
bital, remain qualitatively the same even with the main 
terms of the AL CVC; second, that the resistivity away 


from the AF QCP can be almost well described by taking 
account of only the E CVC. 

Then, I turn to p ab near the AF QCP. From its tem¬ 
perature dependence shown in Fig. |15[ we see three main 
features about p a b in the MT+AL CVC case. First, p a b in 
the MT+AL CVC case shows the T-linear dependence, 
which is similar for the other three cases . This origin 
is the same for the other three case d — i.e. the T 0 5 
dependence of the unrenormalized QP damping of the 
d xz / yz orbital around k = (|^7r, |^7r), since the CVCs 
little affect the power of the temperature dependence of 
Pab and since the d xz + d yz component remains domi¬ 
nant even with the main terms of the AL CVC [see Fig. 
[I6|[a)]. Second, the values of p a b in the MT+AL CVC 
case at high temperatures are nearly the same as those 
in the E CVC case at the corresponding temperatures. 
This is due to the same reason as that away from the AF 
QCP. Third, as temperature decreases, the value of p a b 
in the MT+AL CVC case approaches the value in the 
MT CVC case. This can be understood by combining 
two facts that the MT and the AL CVCs separately be¬ 
comes non-negligible at low temperatures, and that the 
AL CVC near the AF QCP is negligible compared with 
the MT CVC in the presence of the even-parity symmetry 
and rotational symmetry. The mechanism of the former 
fact explained above, and the mechanism of the latter 
explained by the authors of Ref. 37; The explanations 
about the latter are as follows: When the system ap¬ 
proaches an AF QCP characterized by spin fluctuation 
at q = Q qcp, that spin fluctuation gives the leading 
contributions to the MT, AL1, and AL2 CVCs through 
ImV)ff] n (fc — k') in Eq. |7l| and in 


Eqs. (|75j) and ([76]), respectively. Then, although the MT 
CVC becomes more important near the AF QCP, the 
AL1 and AL2 CVCs become little important compared 
with the MT CVC near the AF QCP due to the cancel¬ 
lation between the contributions from k ! and — k ! arising 
from the spin fluctuation at q = Q qcp- This cancella¬ 
tion is because in the terms of the AL1 or AL2 CVC only 
A 1 /; 2 ;cc(fc , ; 0) is odd about momentum (i.e., the others are 
even) due to the even-parity symmetry and because the 
states at — k ' + Qqcp and —fc' — Qqcp are equivalent 
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FIG. 16: Temperature dependence of (a) a xx and orbital-decomposed components in the MT+AL CVC case at U = 2.1 eV 
and the orbital-decomposed components of (b) the d xz and d yz orbitals and (c) d xy orbital in the four cases at U = 2.1 eV. 



T(eV) 


FIG. 17: Temperature dependence of i?H at U = 1.8 eV and 
Jh = §- in the four cases. 


due to the rotational symmetry. 


Moreover, we determine the role of each t 2 g orbital 
in determining p a b near the AF QCP from the results 
of the orbital-decomposed components of a xx . Figure 
16 'a) shows that the main terms of the AL CVC keep 
the dominance of the d xz + d yz component unchanged. 
Furthermore, from Figs. |16[ b) and 16 c), we see a similar 
behavior for p a f, at low temperatures, i.e. an approach of 
the value of the d xz + d yz component or the d xy compo¬ 
nent in the MT+AL CVC case towards that in the MT 
CVC case with decreasing temperature. 


Combining the results at U = 2.1 eV, we find that 
the T-linear p a b and the dominance of the d xz / yz orbital 
which are obtained in the MT CVC case are qualita¬ 
tively unchanged even in the MT+AL CVC case, and 
that there are two almost distinct regions of the temper¬ 
ature dependence of p a b- Those regions consist of high- 
temperature region, governed mainly by the E CVC, and 
low-temperature region, governed mainly by the E CVC 
and the MT CVC. 


2. Hall coefficient 


We start discussing I?h(= o-~ x lim^^o ^jj~) away from 
the AF QCP. We show its temperature dependence in the 
four cases in Fig. and see from that figure three main 
features in the MT+AL CVC case. First, at high tem¬ 
peratures, the values of i?n in the MT+AL CVC case are 
close to the values in the E CVC case. This origin is the 
same for p a b at high temperatures, i.e., the small effects 
of the MT and the AL CVCs due to the large QP damp¬ 
ing. Second, when temperature is low, the value of 
in the MT+AL CVC case becomes almost the same as 
that in the MT CVC case. This result can be understood 
that the main effects of the MT CVC on at low tem¬ 
peratures remain leading even including the main terms 
of the AL CVC. Its mechanism is as follows: As shown 
in Ref. 1.181 the main effects of the MT CVC on i?H at 
low temperatures are the decreases of the negative-sign 
contributions of the d xz + d yz component of the trans¬ 
verse conductivity around k = (|| 7 r, || 7 r) as a result of 
the magnitude changes of the currents of the d xz / yz or¬ 
bital around k = (|^ 7 r, §^ 7 r) due to the MT CVC arising 
from the spin fluctuations of the d xz / yz orbital around 
q = Qic-af- Although the currents of the d xz / yz orbital 
around k = (^ 7 r, are affected by the AL1 and AL2 

CVCs arising from the above spin fluctuations, the main 
effects of the MT CVC persist in the MT+AL CVC case 
due to the cancellation of those AL1 and AL2 CVCs in 
the presence of the even-parity and the rotational symme¬ 
try. (The mechanism of this cancellation 37 ^ was explained 
in Sect. Ill C 1.) It should be noted that we can under¬ 
stand why only for i?n the main effects of the MT CVC 
survive at low temperatures even with the main terms of 
the AL CVC by considering the difference between the 
important factors for <j xx and lim^^o '■ since the im¬ 
portant factor for a xx is the unrenormalized QP damp¬ 
ing, the effects of the MT CVC on the currents of the 
d xz / yz orbital around k = (§ 5 + ^ 7 r) are little important 
for p a b away from the AF QCP due to the large unrenor¬ 
malized QP damping; on the other hand, since not only 
the unrenormalized QP damping but also the momentum 
derivative of the angle of the current becomes important 
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(a) (b) 




FIG. 18: Temperature dependence of limH-s-o and orbital-decomposed components in the MT+AL CVC case at U = 1.8 
eV and the orbital-decomposed components of (b) the d x , and d yz orbitals and (c) d xy orbital in the four cases at U = 1.8 eV. 



T(eV) 


FIG. 19: Temperature dependence of Ah at U = 2.1 eV and 
Jh = jr in the four cases. 


for lim^-yo a ^ y • the effects of the MT CVC on the cur¬ 
rents of the d xz / yz orbital around k = (|^7r, r) become 
considerable for Ah due to the large momentum deriva¬ 
tive. Third, three specific features of Ah obtained in the 
MT CVC case survive even including the main terms of 
the AL CVC; the three specific features are emerging a 
peak, crossing over zero, and increasing monotonically 
in the high-temperature region. This can be understood 
that at high temperatures the effects of the AL CVC are 
small due to the large QP damping, and that the main 
effects of the MT CVC at low temperatures remain qual¬ 
itatively unchanged. 


Next, we analyze the orbital-decomposed components 
of lim ++ away from the AF QCP to determine the role 

ff-s-0 H 


of each t 2 g orbital. Due to the same reason for a xx , we 
consider only the d xz + d yz component and d xy compo¬ 
nent of linitf^o the former and the latter are defined 
as the equations that ]C{o}=i Eq. (39) are replaced 
by, respectively, E{ a }=i and E{ a }= 3 - From Fi S s - 11%)“ 


18’c), we find three main properties: at high tempera¬ 


tures, the effects of all the CVCs on those components 
are very small; in the low-temperature region, the tem¬ 
perature dependence in the MT+AL CVC case is qual¬ 
itatively the same for the MT CVC case; a peak of the 


d X z + d yz component survives even with the main terms 
of the AL CVC. Those results can be understood in the 
similar way for Ah- It should be noted that the difference 
in the important factor is the origin why the d xy orbital 
gives considerable contributions to only lim^^o aF 
though its contributions to a xx are negligible. 


Thus, we deduce from the results at U = 1.8 eV that 
the qualitative behavior of Ah away from the AF QCP 
can be captured by taking account of the £ CVC and the 
MT CVC. 

Then, we go on to analyze the temperature depen¬ 
dence of Ah near the AF QCP. The results at U = 2.1 
eV in the four cases are shown in Fig. [19] We see that 
even including the main terms of the AL CVC, Ah shows 
three specific features (emerging a peak, crossing over 
zero, and increasing monotonically in high-temperature 
region), and that Ah has two almost distinct regions as 
a function of temperature. The former result can be un¬ 
derstood in the same way as the explanations about the 
third result of Ah away from the AF QCP, and the lat¬ 
ter can be understood in the same way for the similar 
property of p a b near the AF QCP (see Sect. Ill C 1). 


Moreover, we determine the orbital dependence of 
liniH->.o Cr fj L near the AF QCP by the analyses of the 
temperature dependence of the d xz + d yz component and 
d xy component. From Figs. [20|(a)-p0|(c), we see three fea¬ 
tures the same as those obtained away from the AF QCP. 
Furthermore, comparing the results away from and near 
the AF QCP, we find that the difference in the d xz + d yz 
component or d xy component between the MT CVC case 
and the MT+AL CVC case becomes smaller near the AF 
QCP than away from the AF QCP. This results from the 
more importance of the MT CVC near the AF QCP. 


Thus, the results at U = 2.1 eV show the validity of the 
qualitative behaviors of Ah near the AF QCP in the MT 
CVC case and the existence of the two almost distinct 
regions of the temperature dependence of Ah near the 
AF QCP, which are the same for p a b near the AF QCP. 
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(a) (b) 




FIG. 20: Temperature dependence of limH->o and the orbital-decomposed components in the MT+AL CVC case at 
U = 2.1 eV and the orbital-decomposed components of (b) the d xz and d yz orbitals and (c) d xy orbital in the four cases at 


U = 2.1 eV. 


IV. DISCUSSIONS 

In this section, I compare the results obtained away 
from or near the AF QCP with several experimental or 
theoretical results. The discussions in Sec. IV A are 
for the comparisons with several experiments, and the 
discussions in Sec. IV B are for the comparisons with 
other theories. 


A. Comparisons with experiments 

In this section, we compare the results obtained away 
from or near the AF QCP with several experiments 
of Si' 2 Ru 04 or Sr 2 Ruo. 975 Tig 025 O 4 , respectively, and 
see that the obtained results are qualitatively consis¬ 
tent with these experiments. In the comparisons with 
Sr 2 Ru 0 . 975 Ti 0 . 02 . 5 O 4 , I believe that the physical origins 
of several behaviors observed experimentally can be de¬ 
duced by comparison with the results obtained near the 
AF QCP. This is because the main effect of the Ti- 
substitution may be the system approaching towards the 
AF QCP compared with S^RuCH; this main effect can 
be treated by increasing the value of U in the model of 
Sr 2 Ru 04 . Although the microscopic mechanism why the 
Ti-substitution causes approaching towards the AF QCP 
is unclear, we can qualitatively understand several dif¬ 
ferences between Sr 2 Ru 04 and SmRuo. 975 ^ 0.02504 as a 
result of the difference in the distance from the AF QCP, 
as I will show below. 

We begin with the comparisons about the magnetic 
properties. The enhancement of the spin susceptibility 
at q = Qic-af away from the AF QCP agrees with the 
neutrorPS or the nuclear-magnetic-resonance (NMR^PS 
measurement of SmRuC^, and the similar enhancement 
near the AF QCP is consistent with the neutron 20 or 
the NMR?° measurement of Sr 2 Ru 0 . 975 ^ 0 . 02504 . Also, 
no sizable commensurate ferromagnetic spin fluctuation 
obtained away from the AF QCP is in agreement with 
the neutron measurement 21 of Sr 2 Ru 04 . Thus, several 
magnetic properties can be well described in the FLEX 
approximation, as explained in Ref. [33J It should be 


noted that to discuss the anisotropy between the inplane 
and the out-of-plane spin susceptibilities, the spin-orbit 
interaction is necessarjRD 

Then, we turn to the comparisons about the electronic 
structure. As discussed in Ref. [35] the larger mass en¬ 
hancement for the d xy orbital than for the d xz / yz orbital 
away from the AF QCP is consistent with an experi¬ 
ment 13 14 in Sr 2 Ru 04 . In addition, the topology of the 
FS away from the AF QCP agrees with the measurement 
of the dHvA effeclP^ or the angle-resolved photoemission 
spectroscopjPsI. However, there is a quantitative differ¬ 
ence in the location of the FS sheet of the d xy orbital near 
k = (| 7 r, | 7 r): in my result, that FS sheet is very close 
to the inner FS sheet of the d xz / yz orbital [see Fig. |9j[a)]; 
in the experiments 13 ^, that FS sheet is not very close 
to the inner FS sheet. That difference exists even in the 
LDA MHZ! 

as described in Ref. [35] To improve that dif¬ 
ference, the spin-orbit interaction of the Ru ions will be 
necessary since the small spin-orbit interaction leads to 
the weak hybridization of the bands near k = (| 7 r, 2 7 r]pS]. 
Actually, that difference is improved in the local-spin- 
density approximation^ including the spin-orbit inter¬ 
action. Although that result indicates the importance of 
the spin-orbit interaction for quantitative discussions, the 
qualitative agreement of the FS is meaningful since that 
qualitative agreement suggests that the present theory 
can capture the aspects of many-body effects on the FS 
of Sr 2 RuC >4 at least on a qualitative level. In addition, 
that qualitative agreement supports the suitability of the 
expectation that the electronic structure of Sr 2 Ru 04 in 
the LDA^ may be regarded as a good starting point 
to include many-body effects. This is another meaningful 
aspect of that qualitative agreement since electron corre¬ 
lation sometimes modifies the FS drastically, resulting in 
the deviation of the FS from the experiment even on a 
qualitative level' 2 . From those comparisons, we deduce 
that the FLEX approximation can qualitatively well de¬ 
scribe the electronic structure of Si^RuCH. 

Finally, I compare the obtained results with the ex¬ 
perimental results about the transport properties. As 
described in Ref. (35] the T 2 dependence of p a t, at low 
temperatures and the non-monotonic temperature de- 
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pendence of i?H away from the AF QCP are qualita¬ 
tively consistent with the experiments- 1 - 1 ^ of p a b and 
i?H of Sr 2 Ru 04 , and the T-linear p ab near the AF QCP 
can qualitatively explain the experimental result 1 ' in 
Sr2Ruo.975Tio.o2504. (The detail of the temperature de¬ 
pendence of i?H hr Sr 2 RuC >4 was described in Sec. I.) Al¬ 
though the spin-orbit interaction generally leads to an ad¬ 
ditional contribution to Ra through the anomalous Hall 
effect-® it has been experiment ally confirmed that such 
contribution is small 19 . Thus, neglecting the spin-orbit 
interaction will be sufficient for at least qualitative dis¬ 
cussions about ■ Combining those discussions, we find 
that the successful descriptions of the transport proper¬ 
ties qualitatively hold in the FLEX approximation with 
the X CVC, the MT CVC, and the main terms of the AL 

cvc. 


B. Comparisons with other theories 

In this section, we compare the results of this paper 
with other theoretical studies and show several better 
points of this theory. First, we focus on the comparisons 
with the DMFT®^ f or a model of Sr 2 RuC> 4 , and show 
that several electronic properties can be better described 
in the method I used due to a sufficient treatment of the 
spatially modulated spin fluctuations. Then, we compare 
my results with the transport properties of Sr 2 Ru 04 ob¬ 
tained in the phenomenological Boltzmann theory within 
the relaxation-time approximation 76 ^, and point out the 
importance of the MT CVC arising from the spin fluctua¬ 
tions of the d xz / yz orbital in order to naturally obtain the 
non-monotonic temperature dependence of Ru without 
any ad hoc parameters of the unrenormalized QP damp¬ 
ing. Finally, we discuss the similarities and differences 
between the main effects of the AL CVC on the charge 
transports in my case and the cas^3 of the single-orbital 
Hubbard model near an AF QCP on a square lattice, and 
conclude that the existence of the two almost distinct re¬ 
gions of the charge transports near an AF QCP as a 
function of temperature is one of the important findings 
of this paper. 

We begin with the comparisons about the magnetic 
properties, the orbital dependence of the mass enhance¬ 
ment, and the modification of t he FS sheets due to 
electron correlation in the DMFT®® for a model of 
Sr 2 Ru 04 in more detail than Ref. ;33, First, in the 
DMFT, spatial correlations, the momentum dependent 
fluctuations, are completely neglected? 7 ', and the effects 
of only local correlations are taken into account 78 ,9 . On 
the other hand, I have shown that not local but spatially 
modulated spin fluctuation with q = Qic-AF is impor¬ 
tant to discuss the magnetic properties of Sr 2 Ru 04 . As 
explained in Sect. IV A, the neutrorP^ and the NMRP^I 
measurement for Sr 2 Ru 04 have observed the enhance¬ 
ment of that spatially modulated spin fluctuation. Thus, 
the magnetic properties can be better described in my 
case. Second, in the DMFT, the locations of the FS 


sheets are more drastically modified due to electron corre¬ 
lation than in my case: on k x = k y line, the 7 sheet of the 
d xy orbital becomes the most inner sheet in the DMFT®*, 
while the most inner sheet in the LDA ® 127 or my case is 
the /? sheet of the d xz / yz orbital. Since the experimental 
result J 19 EH are consistent with the result in the LDA® 2 ?' 
and my case, the agreement about the FS is better in 
my case than in the DMFT®*. Third, in the DMFT, the 
mass enhancement of the d xy orbital is larger than that 
of the d xz / yz orbital for finite values of e.g., the 

former and the latter become 3.2 and 2.4 at U = 2.3 eV 
and J H = 0.3 eV (~ 0.1317) or 4.5 and 3.3 at U = 2.3 eV 
and Jh = 0.4 eV (~ 0.17417). Thus, the larger mass en¬ 
hancement for the d xy orbital than for the d xz / yz orbital 
is obtained in both the DMFT® and the FLEX approxi¬ 
mation, although the large mass enhancement of the d xy 
orbital is realized in a wider region of the parameters 
space in my case than in the DMFT®!; in particular, that 
larger mass enhancement of the d xy orbital is obtained 
even at Jh = 0 eV in my cas^H. it should be noted that 
it is important and necessary to check whether spatial 
correlations, neglected in the DMFT®*, keep the orbital 
dependence of the mass enhancement qualitatively the 
same since the spatial correlations sometimes drastically 
change the results obtained in the DMFT (e.g., see case^E* 
of a single-orbital Hubbard model on a square lattice). 
Actually, even in a two-orbital Hubbard modePl* 0 n a 
square lattice, spatial correlations included by consider¬ 
ing a cluster cause the almost perfect disappearance of 
the unusual Jh dependence of a critical value of U for a 
Mott transition, which is obtained in the DMFT. Com¬ 
bining the discussions of this paragraph, I conclude that 
several electronic properties of Sr 2 RuQ 4 can be better 
described in my case than in the DMFT®®*. In partic¬ 
ular, it is significant to find the importance of the MT 
CVC for obtaining the non-monotonic temperature de¬ 
pendence of 7 ?h of Sr 2 RuC >4 in this paper since in the 
DMFT the CVCs are neglected due to lack of the mo¬ 
mentum dependence of the self-energ d 82 l 83 l. 

I turn to the comparisons about the transport proper¬ 
ties of Sr 2 Ru 04 with the phenomenological Boltzmann 
theory in the relaxation-time approximatiorP^l, in which 
all the CVCs are neglected^. In the relaxation-time ap¬ 
proximation for Sr 2 Ru 04 , the unrenormalized QP damp¬ 
ing, r” 1 , is given by t~ 1 = q a + a a T 2 with the ad- 
hoc parameters r) a and a a , which are chosen as qd xs = 
Vd yz = 2.75, q dxy = 3.25, a d xz = 0.035, a dy , = 0.04, and 
a d = 0.06; this form of Tty 1 as a function of tempera¬ 
ture is typical of the FIl 9 * 10 * 84 *. In addition, only for the 
calculation of the in-plane resistivity, the authors of Ref. 
US added 0.6T to r d 1 since they assumed that Sr 2 RuC >4 
were close to a ferromagnetic instability. Although that 
assumption is experimentally incorrect 21 * 69 !, their results 
at low temperatures will remain qualitatively unchanged 
since the contribution from the d xz / yz orbital is more im¬ 
portant than that from the d xy orbital due to the smaller 
Tty 1 of the d xz / yz orbital. Adopting the phenomenologi¬ 
cal Boltzmann theory in the relaxation-time approxima- 
























26 


tion with those expressions of the unrenormalized QP 
damping to p a b and -Rh, the authors of Ref. [76] obtained 
the T 2 dependence of p ab at low temperatures and the 
non-monotonic temperature dep enden ce of Rn, which are 
consistent with the experiment j 12 l 42 l However, as they 
pointed out in Ref. [76], the result of i?n is very sensitive 
to the small relative variation of and ad z , and the 
sign change of Ah disappears in some cases. In addition, 
and aj yz should be the same due to the tetrago¬ 
nal symmetry of the crystal. Thus, although the results 
obtained in the relaxation-time approximation ' 6 seem to 
be reasonable, the validity of the choice of the ad hoc pa¬ 
rameters is unclear. On the other hand, I have shown in 
Ref. [38] without any ad hoc parameters of the unrenor¬ 
malized QP damping that the MT CVC arising from the 
spin fluctuations of the d xz / yz orbital is essential to ob¬ 
tain the non-monotonic temperature dependence of i?n- 
Furthermore, in this paper, I show that the importance 
of that MT CVC remains unchanged even if we consider 
the main terms of the AL CVC. From those arguments, 
I propose the importance of the MT CVC arising from 
the spin fluctuations of the d xz / yz orbital to understand 
the temperature dependence of I?h of Sr 2 RuC> 4 . 

I close this section with the comparisons about the 
main effects of the AL CVC with the cas ^ 7 of the single¬ 
orbital Hubbard model near an AF QCP on a square 
lattice. In that single-orbital case, the authors of Ref. [37: 
analytically or numerically studied the effects of the AL 
CVC on p ab and i?H near the AF QCP where the spin 
susceptibility at q = (n, tt) was most strongly enhanced. 
Then, their analytic study revealed the cancellation of 
the leading contributions in the AL1 or the AL2 CVC 
near the AF QCP, whose details were explained in Sect. 
Ill C 1, and their numerical study considering only the 
contributions of the MT CVC and the AL CVC from the 
states on the Fermi level revealed the qualitative validity 
of the results obtained without the AL CVC. However, 
due to neglecting the other contributions near the Fermi 
level, they did not obtain the existence of the two almost 
distinct regions of the transport properties near the AF 
QCP as a function of temperature, which is revealed in 
this paper. It should be noted that I obtained the similar 
results®^ to those they had obtained, a decrease of the 
value of the inplane resistivity due to the AL CVC and a 
larger decrease of that at high temperatures than at low 
temperatures (see Fig. 15). Thus, the aspects of the AL 
CVC in my case are qualitatively consistent with those 
in the single-orbital casJ®l In addition, it is one of the 
important findings of this paper to reveal the existence of 
the two almost distinct regions of the charge transports 
near an AF QCP. 


V. SUMMARY 

In summary, after explaining the formal derivations of 
Pab , 7?Hj and the FLEX approximation with the X CVC, 
the MT CVC, and the AL CVC, I studied p a b and Ah 


for a f 2s -orbital Hubbard model in a PM state near or 
away from the AF QCP on a square lattice in the FLEX 
approximation with the X CVC, the MT CVC, and the 
main terms of the AL CVC, and then found the three 
main results about many-body effects. The first main re¬ 
sult is showing that the results of the previous studied ^ 
remain qualitatively unchanged even with the main terms 
of the AL CVC. This indicates the qualitative validity of 
the arguments in the previous studies 10 ^. The second 
main result is finding the two almost distinct regions of 
the charge transports near the AF QCP: p a b and Ah in 
the high-temperature region are described by including 
only the X CVC, while the X CVC and the MT CVC 
are necessary for their descriptions in low-temperature 
region. The third main result is clarifying the differ¬ 
ence of the effects of the MT CVC between p ab and Ah 
away from the AF QCP at low temperatures: the MT 
CVC leads to the considerable effect only on the latter, 
although at high temperatures only the X CVC affects 
p a b and i?H- Thus, the second and third main results 
highlight the important aspects of many-body effects on 
the charge transports. I also showed several results of 
the magnetic properties and the electronic structure for 
that model in the FLEX approximation. Those sup¬ 
p ort a deeper understanding than in the previous stud- 
ie d 10 l 38 l Then, comparing the obtained results with sev¬ 
eral experiments, I achieved the qualitative agreement 
about the momentum dependence of the spin fluctuation 
in Sr 2 Ru 04 21 ^or Sr 2 Ruo. 975 ^ 0 . 025 C> 4 20 ^, the orbital 
dependence of the mass enhancement in Sr 2 Ru 04 1 ®^, 
the topology of the FS of Sr 2 RuOj 131281 , the T 2 depen¬ 
dence 12 of pab at low temperatures in Sr 2 Ru 04 , the non¬ 
monotonic temperature dependence 42 of Ah hr Sr 2 RuC> 4 , 
and the T-linear p ab in Sr 2 Ruo. 975 Ti 0 .o 25 04 1 '. Further¬ 
more, by the comparisons with other theories, I showed 
several stronger points to discuss the electronic proper¬ 
ties of Sr 2 RuC >4 than other theoried^^, and clarified 
the similarities and differences of the main effects of the 
AL CVC between the present nmltiorbital case and a 
single-orbital casd®l 

Several important issues remain for further study. One 
is the extension of the present method to case with the 
spin-orbit interaction. In particular, it is desirable to 
study the anisotropy between the inplane and the out- 
of-plane susceptibility, the quantitative effects on the FS 
deformation including many-body effects, and the charge 
Hall effect including both the usual Hall effect and the 
anomalous Hall effect. It is a lso important to extend the 
present method to cas d 181191 with the Ru0 6 distortions 
near the ferromagnetic QCP since its results and the 
present results lead to deep understanding of the sim¬ 
ilarities and differences between many-body effects on 
the charge transports near the ferromagnetic and the 
AF QCPs. Furthermore, the study about the charge 
transports of the 3D ruthenateiP^ using the extended 
method is useful to clarify the role of the dimensional¬ 
ity in the charge transports of a correlated multiorbital 
system. Then, another important issue is the extended 
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(a) (b) 



FIG . 21: Contours used for the analytic continuations of Eq. 
(101. Contour shown in panel (b) is one of the possibilities. 


study about the charge or heat transports of Sr 2 Ru 04 
in a superconducting phasd^d^ZI s i nce it may clarify the 
role of each i 2ff orbital in the superconducting phase. Fi¬ 
nally, it is challenging and important to study the charge 
transports of other transition-metal oxides 4 or the heavy- 
fermion materials such as CeCoIns^ and UPt or the 
organic conductor:^ by extending the present method 
and adopting that to other nrultiorbital Hubbard models 
or the nrultiorbital Anderson models. Their achievements 
may provide deeper knowledge about ubiquitous or char¬ 
acteristic properties of many-body effects on the charge 
transports of a correlated multiorbital system. 


Appendix A: Derivation of Eq. (13) 


In this appendix, we derive Eq. (131 from Eq. (10). 


In this derivation, we use the analytic properties of the 
single-particle Green’s function and reducible four-point 
vertex function in terms of frequency variables. The for¬ 
mer becomes singular for line = 0 with e being its fre¬ 
quency variable. The analytic property of the latter is 
the same for the two-particle Green’s function^ it be¬ 
comes singular when the frequency variables e, e', and 
ui satisfy line = 0 or (Ime + Imw) = 0 or Ime' = 0 
or (Ime' + Imu) = 0 or (Ime + Ime' + Imw) = 0 or 
(Ime — Ime') = 0, where we consider case for Imu > 0 to 
take ifl n + i0+. As a result, there are 16 possibil¬ 


ities of the four-point vertex function in real-frequency 
representation, as shown in Table [i] 

Using those analytic properties, we can carry out the 
analytic continuations of the first and the second terms 
of Eq. ( |10[ ) . This is because we can rewrite the summa¬ 
tion about the Matsubara frequency as the correspond¬ 
ing contour integral®. Using the contour G shown in Fig. 
21 ’a), we can carry out the analytic continuation for the 


first term: 
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where —> represents the replacement of i£l n by u> + *0+, 
and gi-acdb{k\u) are 


gv,acdb(k-,co) = G™(k,e + u>)G™(k,e), (A2) 

g2;acdb{k;uj) = G$ (k, e + w)G^ } (fc, e), (A3) 

and 

93;acdb{k-,u) = G^ (k, e + co)G^\k, e). (A4) 

Next, by replacing TJ2 m an d in the sec °nd term 

of Eq. ( flo]) by the contour integrals with the contour 
G in Fig. |21[a) and the contour C such as Fig. H3 b ): 
respectively, we can similarly carry out the analytic con¬ 
tinuation of the second term: 
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TABLE I: Relation between the additional subscripts of the four-point vertex function and the inequalities of its frequency 
variables. Because of Imw > 0 and the analytic properties of the four-point vertex function, there are 16 possibilities. 
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^ ^ ^ 'j g'ki-’^ba i'^k'v)cd, ^ ^ 


de 


k,k‘ {a} 




-oo 47r * 7-00 47r * 


de' e' 

— tanh — 

2 T 


x ] (tanh^ -tanh^^)(/ 2 ;aAB&(fc;w) r 21;{A} (ft, fc'; u)gi ;Cc dD{k';ui) - r 22 . II;{A} (fc, k';u)g 2 -,ccdD{k!-, w) 

e + w 


{( ‘ e+ , 

- tanh -g 3 -aABb(k ; w) r 31 _ I;{A} (fc, u)gi-CcdD(k' ; w) - r 32;{A} (fc, k'-,u)g 2 .ccdD(k'; w) 
+ tanh ~9i-aABb{k-, w) r n _ I;{A} (fc, k';uj)g 1 . Cc dD{k'\ u) - T 12;{A} (ft, k!\u)g 2 . C cdD (ft'; w 


1 

TV 2 


'^2'^2{Vkv)ba{Vk'v)cd'^2, [ 47ri [ 

fe.fe' {a} {A} 1 ^ 00 


/ °° E ta nh e ' + W 

— OO 


47ri 


2T 


x <j (tanh ^77 - tanh g 2;a ABb(k\ ui) T 2 2-iv-,{ A} {k,k';u})g 2 -CcdD(k'-,uj) -T 23 . {A }(k,k'-,u)g 3 -c c dD(k';u}] 


{(t 2T 

- tanh —^-g 3l aABb{k', w) r 32;{A} (fc, k';u})g 2 -CcdD{k'; u) - r 33 _ I;{A} (fc, k';cj)g 3;C cdD(k'; u) 
+ tanh ^ 9i-aABb(k;cd ) ri 2;{A} (fc, ft'; uj)g 2 -CcdD{k'] w) - r 13 _ I;{A} (fc, fc'; uj)g 3; CcdD{k'; w 


N 2 


^ ^ 'fa kv^baiyk'v)cd ^ ( 


de Z 00 de' , e' + e + w 
coth 


{a} 


{^4} 


-oo 47 ™ 7-oo 4^ 


2T 


x 1 (tanh ^77 -tanh^)5f 2;aAB& (A:;w)3 2 ;C c dn(fc , ;w) r 22 . ra . {A} (fc, fe'; w) - r 22 . IV; {A}(fc, w) 


2 T 
e + uj 


{(‘ 

- tanh ' L ^^g 3 -aABb{k\ w) (r 31 _ II; { A} (fc, k'; w) - r 3 i_i ;{A} (fc, k 1 ; w)J gi-CcdD(k w) 
+ tanh A; gi-aABbik ; u)g 3 . Cc dD{k'\ w) T i 3 _i ; {a} (*, fc'; w) - T i 3 . n;{A} (fc, fc'; w)J | 


de 


de' 


= _ 7V 2 E / ^EEW^(^)“> tanh w) j j—'^2Ju-{A}{k,k'-,uj)gi. i ccdD{k';uj) 


k,k’ 


, „ 47T'j 

{a} {A} 


2T 


-oo ^ /=1 


( 


tanh - tanh g 2 - a ABb(k; w) 


r°° de' 3 


47ri 


T E ^2i;{A} (ft, ft'; u)9hCcdD{k '; w) 


- tanh^-^g 3;aAB& (fc;w) / E E ^ ; {A} (ft, ft'; u)gi-,CcdD(k'; ui) 


de 


2 T 


47ri 


(A5) 
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where Jw^{A}{k 1 k'; ui) are 


£ c — 6 r ■ 

Ju-,{A} (k, k';u>) = tanh — r n . I;{A} (fc, fc'; u) + coth |Tn_ii;{A} (&, w) - r n _i ;{j4} (fe, fc'; w) 


^i 2 ;{A}(fc,fc';w) = (tanh 6 + ^ - tanh ^) r i 2 ;{A}(&, *';<*>), 


(A6) 

(A7) 


e + e' + w 


Ji 3 ;{A}(fc,fc , ;w) = - tanh 2r r i 3 -i;{A} (k, k'\u)- coth 
J2i-,{A}(k,k';uj) = tanh ^r 21;{A} (fc, fc'jw), 


ri 3 -H;{A}(fc, fc';w) - r 13 _ I;{j 4}(fc,/c';w) , (A8) 

(A9) 


J22-,{A}{k,k'\u}) = (coth 


-tanh — )T 22 -n-,{A}(k,k';u)) + (coth -—— - coth — J T 2 2-iii ; {^.} (A , k 1 ; w) 


2T 


+ (tanh 


e' + uj 
2 T 


2 T 
— coth 


e + e + u> 


2 T 

^ | Cc^ 

J 2 z-,{A}{k,k'\u) = - tanh r 23;{A} (fc,fc'; u;), 


r22-IV;{A}(fc, fc 7 ; w), 


e + e' + w 


J3i-,{A}(k,k';u>) = tanh—r 31 _ I;{A} (fc,fc';w) + coth 
J32-,{A}(k,k'-,u) = (tanh 6 + ^ - tanh ^)r 32; { A }(fc, k'\ u), 


r 3 i-n ■{A}{k,k l -,ui) - T 31 . 1 -{A}(k,k'-,uj) 


(A10) 

(All) 

(A12) 

(A13) 


and 


e'+w, 


e — e 


J 33 ;{ J 4 }(fc,fc , ;w) =- tanh -^^r 33 _ I;{A} (fc, fc'; w) -coth ^ 


r 3 3-n ; {A}(fc, *';w) - r 33 _i ;{A} (A;,fc , ;w) 


(A14) 


In Eq. (A5), I have not explicitly written whether the After some algebra for Eq. (Bl), we obtain the Bethe- 


integral is the principal integral (containing a hyperbolic 
cotangent of frequency). 


(141. 


Combining Eqs. (|Al| and (A5), we finally obtain Eq. 


Salpeter equation to determine A u -AB(k, ie m ; q , if2„), 


Appendix B: Derivation of Eq. © 


In this appendix, we derive Eq. (16) after defin¬ 


ing the three-point vector vertex function in Matsubara- 
frequency representation and carrying out its analytic 
continuation. 


We can express Eq. (131 in a more compact form 


since JwAA}{k, k'\uS) are connected with the three-point 
vector vertex function in real-frequency representation. 
We begin with the definitiorJ® in Matsubara-frequency 
representation, A !/;AB (fe, ie m ; q, ifl n ) = A„ ;AB (fc + 
q^ k, idrn) • 

^ G a A{k A q, ic rn -\-n)^i';AB{k^ miQi ift n )G B b(k,ie 

ra) 


A,B 

= [ dre iern+nT I dr'e 

Jo Jo 


( r ^rCk+qa{ T )J-qv( T )c fcb ). 

(Bl) 


J^u\ab{J^i^m'i Qi^^n) — kv)ab 

J 1 

+ AfEim I" abCD^k, it m ? k'itm’] q, ikln) 
k' m' {A} 

X GcA(k -f q, itm'+r^GBDik i ) (Vk f v')AB • (B2) 


Since the second term of Eq. (B2) has the similar ana¬ 


lytic property for the second term of Eq. (10), we can 


similarly carry out the analytic continuation of the for¬ 
mer as follows: 
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(a) (b) (c) 



FIG. 22: Contours used for the analytic continuations of the three-point vector vertex function. line > 0 and line + Imu > 0 
are satisfied in panel (a), Ime < 0 and line + Imcu > 0 are satisfied in panel (b), and line < 0 and Ime + Imu; < 0 are satisfied 
in panel (c). 


1 

itEEE r abCD ^ ^m' i Q1 ^n)^CA(^ Q1 ^m'-\-n)^BD 1 ^m ') i^k'v) AB 


N 

1 

N 


k' m' {A} 


r f de' e' 

EE / . tanh r abCD{k, zc m , A: , e ■ q ■ i^l n ^G qA^ k -f- e/. e z^l n ^G bb (A: , e ) 

t' 2 C" 


fc' {A} 


d - (IT abCD{k^ ^e m , k , te m , q, i£L n )GcA{k T q, -t- i^l n ^)G bd( k , fe m ) 

+ ( k , ze m , k , z^l n . q, zQ^G c a{ k -t- 17 , itm)G bd ; te m 2 fl n ) {y^ ^ ab? 


(B3) 


where the contour C" is the contour shown in Fig. 22 'a) 
or [22j' b) or |22]c), depending on the values of Ime and 
Ime + Imw for the second term of Eq. ( |B2[ ). Thus, the 
three-point vector vertex function in real-frequency rep¬ 
resentation is given by 

A-v;l;ab(k : — ( ^ku)ab 

+ Ti E E E [ %-^‘^CDik, k’-q) 


k' {A} l' — l J 

X gi';CABD(k';q)(vk'is)AB , 


(B4) 


where case for l = 1 or 2 or 3 corresponds to case of Fig. 
[22^ a) or |22]b) or 22 c), respectively. Using Eqs. (131, 


114p, and (1IMI), we can rewrite Eq. (13) as Eq. (16). 


Appendix C: Derivation of Eq. (36) 


linear order of q as 

GaA^k —, ic m )G Bb{k-\-i 
^ GaA^k , idm)GBb(k , t£ m _|_ n ) 


+ Ef [G'aA(fe,*e m )^G B 6 (fe,ie m+n ) , (Cl) 






where we introduce a quantity, 


, xZ 1 _ 5/i(x) dh(x) 

g (x)- si h ( x )\=9(*)- Si -(C2) 


Second, A. y; d c (fe±, ie m +„, fc±, «e m ) is approximated as 
^y^dciki: > t€ m _|_ n , fc-j- , fe m ) ~ A y ;dc(Ab 26 m _|_ n , Ag 2C m ) 


i ^ ^ ^y,d,c{k: te m _(_ n , A:, ie m ). 


2 dA, 


(C3) 


In this appendix, I explain the detail of the deriva¬ 
tion of the q-linear terms of K x vl/ (q,itl n ) in the most- 
divergent-term approximation^-®!]. 

We start to show the leading-order q dependence 
of each quantity appearing in Eq. ( [35] ) . First, 
G a A{k-, ie m )GBb(k + , ie m + n ) is approximated within the 


Third, the leading-order q dependence of 
K-AB{k_,ie m ,k + ,ie m+n ) is given by 

A-x;AB (An~,ie m , A;_|_,fe m p n ) ~ A x;AB(k,ic mi ^1 ^ra+n) 

H - ^A-x;AB (^1 i^mi ^1'i^m+n') 1 (C4) 
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with the equation derived from Eq. (B2), 
AA x;j 4 B (fc, ie 

mi ^i ^m+n) 


_ T ^ ^ -pi 1 ) (u ■ u ■ ./ ■ »/ • \ since we can negiecr me q-nnear rerm arising 

~jV ABcdtK>*em,fc,*em+n ; K , *e m <, fc ,*e m '+nj from a pair of tw0 s i n gl e -particle Green’s func- 


vertex functiort- 47 * 48 !. Fifth, the q dependence of 

f g Gfb(k— , iCiji') A.j^-g f (Aj-|-, ie m , k— i ^rn)G C g{k-\-^ ie m ) Is 

negligible in the most-divergent-term approximatioif 47 48 
since we can neglect the g-linear term arising 


k' m' {a} 


]T y G ca {k' ,ie ml )-^-G bd {k' ,ie m ' +n ) 


x > ^ 


2 L 






X ^-x\abi}z j 5 ^ 5 ^m'+n) 

T ^nt^-v^(i) 


tions whose frequencies are the same in the most- 
divergent-term approximation 47 48 . More precisely, 

f g Gfb(k— , ie m ) A^.^y (^+7 7 ^m)G C g(k _^, ie m ) is 

approximated as 


/y 5Z5I5Z r ABcd( fe ’ ie m,k : ie m+n ; fc', ie m /, fc', ie m 7 + „) H Gf b (k-,ie m )A a , gf (k + ,ie m , fc_, ze m )G cg (fc+, *e m ) 

m' {a} 


/>0 


x Gbd(k , ze m /_|_ n )(jr ca (&', ie 

m ') 

x AA x . ab (k ^icm/^k , ie m '_|_ n ). 


(C5) 


In Eq. (C5), we have neglected the q-linear term 
arising from the irreducible four-point vertex function 
because that is negligible compared with the main 
terms in the most-divergent-term approximation- 7 ®. 
Fourth, due to the same reason, we can neglect the 
q-linear term arising from the irreducible six-point 


E Gfb(k, ie m )A a - tg f(k, ie mi ki ^m)^cg (^7 ^ m) 
f i9 

dG c b(k, ie 

m) 

9/c„ 


(C6) 


In the final line, we use one of the Ward identitie d - 47 -! — l 
Using the q dependences of the quantities appearing 
in Eq. (35), we obtain the q-linear terms of K xvu {q, ifl n ) 
in the most-divergent-term approximation 4 ^® 
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jn ^ ^ ^ 

K xyv {q,iCl n ) = ‘WjyEE E Qk 

k m a,b,A,B * 


d(vky) 


ba 


GaA^k, (^? ^? ^m+n)^B6(^j ^m+n) 


T 


d( v ky) 


6 ^y N J2Yl E dk 

k m a.b.A.B 


T 


^m+Ti) ^ ^ GsbiJ^i ^ 771 +n) 




^tEEE aa x\ba{}^i it m? k"> ^m+n)G ad(k ? ^ m+n) 


T 


k m {a} 


2 L 

v 

dG c b{k, ie 

m) 

dk u 


^-y',dc{k") ^m+n? ^? ^m) 


^EEEa x;ba(k , JCm) *£m+n) ^ ' 


Qrj ^ad(^i '^m+n) ^G c fc(A;,ie m ) 


T 

TV 


fc m { a } 






Ay;rf c (fc, Te m _|_ n , k : idrn) 


EEE A * ;ba(fc,*e m , k , ze m _|_ n )G^ ac ^(fe, ze m _|_ n ) 


^G c fr(/c, Z€ m ) V - ^ Qy d^-y,dc{k-> ^m-\-ni k, i^-rn) 


k m { a } 


<9£+ 


E 


<9fc„ 


^EEE AA x . ba (k , Te m , fc, ie m+n )G cb (k , ie m ) ^ Gad ^d €m+n ') A y . dc {k , ie m+ni k, ie m ) 

fc m {a} ^ 


iy ^ ^ ^ ^ ^ ^ h-x-,ba(k-> itn 2 ? &? ^ e m+n) ^ ^ 


Qrj dGcb{.k)i€m) ^^ad (^5 ^m+n) 


AT 


T 

-4EEE A x;ba (k, 26 m 5 k , ^Gm+n)G c fr (fc, ^ m) ” 


fc m {a} 


dk„ 


dk„ 


Ay-dc{k , Te m _|_ n , k : iCrn') 


dG ad (k, Tc m -(- n ) v ' dAy^dc^k, iCm+ni k 7 iCm') 


N 

k m { a } 

' T \ 3 




E 


<9fc„ 


(*) E E EEE Gsb{k , ie m )AA x . ba (k , *e m , k , ie m+n )G a ^(fc, ze m + n ) 

k,k',k" m,m',m" {a} {A} -F,G 

dG GF (k n 


0*:: 


G Dd(k , Z6 m / _|_ n )Ay|^ c (fc I'itrr I'+ni ^ ; ^m / )^’cc(^ ?^m / ) 


X ^3 ; ABC-DFg(^’ m+n? kj ^m] & ? ^m'? ^ -\-n->k ^iCui'f^k ,ZC m ") 

^_x 

rp g ' (j 

-(+) E E EEE A-x-,ba(ki is m 1 k , Z£m+n )E f [G Sb (fc, ie m ) — G aA (k, ie m+n ) 

k,k' ,k" m,m' ,m" {a} {A}-F,G rj ^ 

dG G F{k" ,ie m ") 


dk" 


G Dd,{k , Te m ' + n ')Ay f ic{k , Te m '-|_ n , A; : i^m')G c c{k , Te m ') 


X -^3;ASCDFG(A ; '*e m+n? ^? 'itm’i k , Te m ', k , Te m '_|_ n , k , T6 m ^, k , T6 m ^). 


(C7) 


The above ninth and tenth terms can be rewritten in a relation® between the irreducible four-point and the ir- 


simpler form by using Eqs. (B2) and (C5), the exchange reducible six-point vertex function, 
symmetry 48 of the four-point vertex function, and the 


/ d d \ 

\ K ~0k~ A J Te m , , Te m '_f_ n , k , Te m ') 


9fc; 

T 


= irEEE ^3;ABCDFg(^’ m+n? k , ^ m\k , it-nri' ? k , i£ m '_|_ n , /c , Z£ m n , , Z£ m " ) 


fc" m" F,G 


9G G F(fe // ?^n 

dk" 


(C8) 


Namely, the ninth and tenth terms become 
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(*) EE EE G Bb (fc, ie m ) AA x;ba (k, ie mi k, ^m+n)^aA(^j ^ m+n )G Dd{}^ i ^m'+n) 

m,m/ {a} {A} 

X Ay-d c (k , Z6 m /_|_ ? 7,, k , 26777' )^cC 5 ^m ') ^ “ 1 “ ^ ^ AS DC (^’ k , 26777, A , Z6 m /_)_ n , A , 26777' ) 


w) EEEEa x\ba(ki it mi k , ze m -|_ n )E f[G Bb (fc,*e 

m )-r^-G aA {k,i^ m+n ) | GDd{k , ZCTTiZ-i-y},) 

fc,fc' m,m' {a} {A} 

(9 5 


X A y; dc(k , 26777' _| -mk , ZEm'^G c C (A “1“ j ^'aBDC ^ 5 ^5 ^ 1 i^m'-\-ni k ,26777,') 

= ~^7 ^ ^ ^ ^ ^ ^ GBb(ki i^m)^A x fi a (ki lE-mi ki i€m-\-n)GaA(k , i'i'^m+n) T^jC~ ^j/;Ab(^i ^m+ri) ^ ^m) {^ky^AB 

k m a,b,A,B V 

/ 'J 1 \2 

-+j EEEE G Bb (fc,ie m )AA x;ba (fc,ie 777, A, 26777-|_77,)G^ aj 4 (fc, 26777,-1-77,) 

fc.fc' m,m/ {a} {A} 

X 9 k f ^*Dd(fe 1 ^m'-\-n)Ay \dc{k , ^6777'_|_77,, k , 'i'€m , ')(*cc(,k , 26777' )J (fc, 26777_)_77, 267775/c ?^m') 

T ^ 

~l~ ^ ^ ^ ^ a;;ba (A^, 26777,/c, 26777-1-77) GBbiJ^i ^6777) G a A^k^ iCm-\-n) 


k m a,b,A,B 

d r, 


X £^A ^2/iAS (A^, 26777_|_77, A, 26777) ky^AB 

- T \ 2 ' 


V 


(+) EEEEa x-ba( k i ie mi A, iEm+fi )E -y- |^G Bb (fc, ie 

m 

^ ( 9 A m+n) 

fc,fc' m,m' {a} {A} V ^ 

d [ 

X GDd{k , l€m'-\-n')Ay ]dc(k , 26 t 77'_)_77 , A , ^m'^)G c c{k , 26777')J A, 26777_)_ 77 , A, 267775 A , 26 m'+n?^ ?^m'] 

T 

N 


k m a.b.A.B 


^ ^ ^ ^ ■j^^A ;E + a (A, ^6777, A, 26777_|_77)£r£?5(A, 'i^-rn)G a A{ki ^m+n) 

' \S 

Ax\ba{J^i i^mi 1 ^m+n) ^ ^ 2 ^Bbi}*i ^m) ^“^aA (^5 ^m+n)j ^ 


k m c,d,C,D 


Ay^ABi,^) i€m-\-ni J^i i€-m) i^ky)AB 

d r 
<9A' L 


SA*, . 

T 9 r 

-*EE E AAa; ;B (^(A; , ie m /_,ze m ^) 77/ G Bc ^(/»! , r i'£m, , -\-n')^-y\dc{k ,ze m /_|_7T,, k 1 xcm'}G c c(k ,ie m 7 ^ 


Thus, replacing the ninth and tenth terms of Eq. (C7) simpler form: 


by the terms of Eq. (C9), we can rewrite Eq. (C7) in a 


(C9) 
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jr _^ ^ ^ ^ 

Kxyi'(Q-) = 1 \t E E. E A^+fc, ie TO , fc, ie m+n ) e 


Qrf dGad(k 5 ^m+n) dG cb {k, idm) 


N 


k m { a } 

T 
N 


dk n 


dk v 


A y -dc{k , ie m+n 5 k, itm) 


EEEa :r;&a {k^it mi k, ie m _|_ n )G^ a d(fe, ze m _|_ n ) 


&Gcbik-fitiri) Qrj dAy dc {k^ Z6 m _j_ n , k, 2€ m ) 


k m {a} 


<9£+ 


E 


Y y ' y ' y ' A x .;, a (fe, J£ m , k, ie m + n ) ^ ' 


Qrf dG cb (k, iCm) dG ac [{k 1 ic 7n -\- n ') 


N 


k m { a } 


dk„ 


dk u 


dk„ 


A-y;dc(k, it rn j rni Al, 26 m ) 


T \ ' \ ' \ ' . ,, . /, . > 9G a( i{kj ic rn j rn ) \ ' (] /i dA.y \dc{ki itm-\-ni k , ^Cm) 

_ jY / v / v / v ^-x\ba\^i 'Itmi &i 'ltm-\-n)t*cb\Ki ^mj 7^ / ^ 

fc m { a } v ri 77 


^EEEa x;ba(fe,*e mi k, £Cm+n )E y [G c &(Me m ) Gad^k^it m+n) 


T 
N 

k m {a} 

^Ay^ c (fc, ie m _(_ n ,/c, Z6 m ) 




T 


= + EEEa x;ba{k, ie mi k , Z 6 m _(_ n )Ay.^ c (fc, ie m .|_ n , k , ic m ) 


k m { a } 

Qr) ( dGad(kiit-m-\-n) dG c b{kiitm) dG a d(kiitm-\-ri) dG c b{kiitm) 


dkn 


>Ef{ 

^EEEa X ; 6a (fc,j£ m ,fc, 


9^ 


^m+r; 


dk u 


dk„ 


k m { a } 




E|{ G cb (k, ie m )-^^G ad (k, ie m+n ) 


dk v 

E 


dky-dc(k i it-m-\-ni k) it-rri) 


dk v 


Gcb{k, ie m ) ^ Gad{ki it-m-\-n)\ 


9 Ay.dc{k, it-YYi+ni k , iCm) 


<9£+ 


2 y 

= ^ (.Qx^v,y Q.y^u,x) ^ ^ ^ ^ ^ ^ ^-x\ba{ki it-mi k , ze m _|_ n ) 


fc m {a} 




X ^ i^-m) Gad{ki itm-\-n) 




Gcbikiit-m) 7 ^ad(^5 7 ^m+n) ^ 
a/c y 


dA.ydc(yki it-m+ni k-, it-m) 

dky 

dA y . dc (k,ie m-\-m k, it-m) 


dk T 


y [0_x&OL,y Q_y$a,x) ^ ^ ^ ^ ' A x -ba{k , ifrn, k 7 Z6 m _(_ n ) 

km {o} 

/ dG cb {k, 26 m ) dG a( i{k, 'ic m _|_ n ) dG cb [k 1 icm'j dGadik , i£m+?i) 




9fc r 


I 




A-y;dc(k, iCm+n, k, iCm'). (CIO) 


In the above derivation, we have used the fact that the 
surface terms arising from the partial integrations about 
fcj; or A;^ become zero due to the periodicity of the Bril- 
louin zone, while I have not used the replacement used 
in Refs. BTIandBBI 


Finally, we can rewrite Eq. (CIOI as Eq. (361 by using 
the equivalence between the i and the y directions. 


Appendix D: Technical details about the numerical 
calculations of the FLEX approximation 


In this appendix, I remark on several techniques about 
the numerical calculations to solve a set of the equations 
of the FLEX approximation self-consistently by iteration 
using the fast Fourier transformation (FFT+l. 

To use the FFT for the quantities as a function of 
fermionic Matsubara frequency, we need to use the zero 
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padding^. This is because of the antiperiodicity in terms 
of fermionic Matsubara frequency®®. For example, when 
the number of Matsubara frequencies is 2 M, the nonin¬ 
teracting single-particle Green’s function should satisfy 

f G° b (fc,ie m ) for 0 < m < M 
G° ab {k, ie m ) = < 0 for M < m < 3M 

[ G° b (fc,ie m _ 4 M) for 3 M <m< 4M 

(Dl) 

The similar property is satisfied for G a b(k,ie m ) and 
H a b(k, iem). Furthermore, due to this property, the quan¬ 
tities as a function of bosonic Matsubara frequency such 
as Xabcd{q,iQn) satisfy the following property: 


Xabcd(q,ift n ) for 0 < n < 2 M 
Xabcdiq , ifin) = { o for n = 2M 

Xabcd(q,itt n -4M) for 2 M <m< 4 M 

(D2) 


Using the above properties, we can utilize the FFT 
to solve the set of the self-consistent equations of the 
FLEX approximation by iteration. First, using the in¬ 
put of the self-energy, which is zero for the first iteration, 
we determine G a b(k,ie m ) from the Dyson equation [i.e., 
Eq. (48)]. Second, we calculate the chemical potential 
from Eq. (52) by the bisection method; in this calcula¬ 


tion, to reduce the numerical error arising from the cut¬ 
off frequency, we use Eq. (521 instead of the equation 


where only G a b(fe,ie m ) appears, and we set the chemical 
potentials in f(e a (k)), G ab {k,ie m ), and G° ab (k,ie m ) the 
same. Third, we carry out the Fourier transformations 
of G a b{k,ie m ) about momentum and frequency: 


1 


4 M -1 


G ab {r 1 T l )=^Y, eikre ^ T E e~ 2 m ^G ab (k,ie m ). 

k 


m —0 


(D3) 


Fourth, using Eqs. (D3) and (441, we determine 


Xabcd(q,i^n) from the equation, 


Xabcd(q > *U n ) — 'y ] 


4M-1 

e - lq -r_}_ XT' e 2 ™™ 


4 MT 


E 

1=0 


x G ac {r,Ti)G db (-r,-Ti). (D4) 


Fifth, from Eqs. (D4), (|45|), and (46), we calculate 


Xn.h,d (q,^n) and Xabcdiq^ri)- Sixth, by solving Eq. 
(51), we obtain V a b c d(q,i^n)- Seventh, we carry out the 


Fourier transformations of Vabcd{q,i^n) about momen¬ 
tum and frequency to determine T, ac (k 1 ie m ) from the 


following equation obtained from Eq. (50): 


4 M —1 

r 1=0 

x ^2V a bcd(r,Ti)G bd (r,Ti). (D5) 

b,d 

At this stage, we obtain the output of the self-energy. 
If the sum of the difference between the absolute values 
of the output and the input becomes less than 10~ 4 , I 
assume that the solution is obtained; otherwise, we re¬ 
place the input of the self-energy by the average of the 
input and the output, and solve the above procedures 
again. Note that if we use the relations about the sym¬ 
metry of the system such as time-reversal symmetry and 
even-parity symmetry and utilize the arrays efficiently, 
we can reduce the memory of the arrays and the time of 
the numerical calculations. 

In the above iterative procedures, we should increase 
the Hubbard interaction terms so slowly as to keep 
the susceptibilities finiteP^ since the calculations are re¬ 
stricted to a PM phase. For example, if we analyze case 
at U = 1.8 eV, J H = J' = and U' = U - 2 J H = 
we begin with U = 0.2 eV and then increase the value of 
U slowly after several iterations. 


Appendix E: Derivation of Eqs. (54)—(57) 


In this appendix, we derive the set of Eqs. (541 (57) 
by adopting Eq. (53) to Eq. (50). 

Using Eq. 


in the FLEX approximation, we ob¬ 
tain T^ bcd (k,ie m ,k'q,iCl n ) in the FLEX approxi¬ 
mation. Setting q = 0 and Cl n = 0 in Eq. (53) and 
substituting Eq. (|50|) into Eq. (|53|), we have 


r ilL(fc, frm, k',ie m r, 0,0) 


= XT E E Va BbD{q' 


N 


q' ,n' B,D 


SG B p{k - g\ie m—n ') 
5G cd (k',ie m >) 


T ST' &VaBbD(q' G^n') n ,, / . n 

2^ E - q , ie m - n '). 


N ^ 6G cd (k',ie m >) 

q ,n' B,D 


(El) 


The first term of Eq. (El) gives the MT term, 

r abed (ki fomi k , 0 , 0) — V a cbd{k k , frm'); 

(E2) 


and the second term gives the AL1 and the AL2 term, 
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r il ) cd L ( fc > ie mi k', ie m r,0, 0) 


abed 

E E E [2 E U aBa'b' 


3 ^ —S fiXa'b’c’d'iQ'’ i^n') TT S 


N 


q' n' B,D {a/} 


SG cd (k',ie m >) 


rrs , 1 ttC ^Xa'b'c'd' (g',*^n') TT c 
u c'd'bD + ~A 2_^t U aBa'b' ^ VW ■ : u ' 


E j-rfl SXa'b'c'd'WG^ri) yrfl 
U aa'Bh> f . ■ N U n> 


{a'} 

l w vr 3 

N 


6G cd {k',ie m ') 


c'bd'D 


{<*'} 

Ggfl(fe - q',*e m—n ') 


SG c d(k',iem /) 


c'd'bD 


H E H U H "•jm'b' ( 9 ' • i!i »' ) w&d.jd («', 

q' n' B,D {A'} C ' ’ m ' 

, 1 A~rC ,„/ „-o ^XA'B'C'D'(q',*n n /) ~ rC „. 0 ■, 

+ 2 2^ N aBA'B'(Q ,*«»') ^ ffe , ■ a N c , D , bD {q >»*V) 


{A 1 } 

-E^ 

{A'} 

T 
~ N 


ti 

aA'BB 1 


SXA'B'C'D' {q ', *^n') r? -t4. 

G < 


5G c d(k',ie m ') 


C'bD'D 


Ggjj^k Q ,2C m _ n ') 


EEE[^E S BcA(g , ,^n')^!c f) i 3 (g , ,^n') + ^a C BcA(9 , ,^')E C C6i ? (9 , ) ^') - 

q' ri {A} 


X GcA(k Q i^m'—ri'jGBD^k Q ,ie m _ n /) 

wvr 3 

N 


EEE [2'^o^Ad(9 , )*^ri')-^C'c6i3(9 , )*^n') + ^aBAdW^ri')N§ cbD {q', i^ri) - Ul\ Bd U ( 

q’ ri {A} 

x G AC (k! + q' ,iem'+ri)G B D(k - q', ie m -ri) 

T 

=-^EEE WaBcA-dCbD^q’^^ri\q'G^ri)Gc A {k' - q', ie m >- n ')G B D(k - q',ie m _ n /) 


■n 

CbcD 


N 

T 


q 1 ri {A} 


-^EEE ^aBAd;C , cfa-D( < 3 ,, > </; ri)G A c{k' + q’, l£ m'-\-ri )G BD {k - q',ie m -ri) 


q' ri {A} 


= r ib ) c d L1 (*. *w; 0 , 0) + r^ L2 (fc, ie m , fc', *e m /; 0, 0). 


(E3) 


In the derivation of Eq. (E3), we have used several equa- and 


tions: the second line of Eq. (E3| is obtained by using 


rrS ^Xa'b'c'd' Wi^ri) TT s 
O' nUn'h' vw : ; o. 


c'd'bD 


5G cd (k’ ,ie m ') 

- V U s lA/Sr 1 (a' m ,\ d XA' B 'C'D'(q ',»»>.') 

- )a'b'A'B'[q^ri) SGcdik ^ iCml) 

X ^ c'D'c’d' (q',in n -)U* 

d'bD 


vs / / -n Axa'B'C'D'W, i^ri) ,%s 


rrC ^ Xg'b'c'd' (q'l i^ri ) rr C 

QBa ' b ' 6G cd (k',iem') c ' d ' bD 
- ST TjC c M C r l / / 0 q ^XA'B'C'J'(g',^n') 

x ^C'D'c'd'iq 1 \i£lri)u3 d'bD 


VC // -o AXA' B 'C'D’{q',i&ri) jr T c 


— N a BA' B '(q' l^ri)- 


SG cd (k\ie 

m ') 


Nc'D'bD^q 'i 


(E5) 


where N^ bcd (q', iCl n ') and N° bcd (q ', are 


Nabcdiq' — ^a,A^b, B + ^2 ^abC DXC DAbW i *^n') , 

C,D 

(E6) 

and 


— NaBA'B'W 


8G cd {k’,i£n 


A - rS . , . Ng bc d{q' , ittri) — 8 a}A S b , B — ^ Uqb CD XcDA B (q'G^ri), 

N C 'D'bD\q > C,D 

(E4) ( E7 ) 
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respectively, and (M s ) a^) 1 c(^ (q , , and 

m C )aL(^>^nO are the inverse matrices of 
M abcdW> iil n>) and M abcd(<l' respectively, 

with 

abcdi 9 = &a,c&b,d ~~ ^ ' XabABjq t ^n')U ABcd , 

A,B 

(E8) 


and 


Mabcdi . 1 1 &a,c&b,d + ^ ' XabABjQ i *Hn')^ABc(ii 

A,B 

(E9) 


to obtain the third line of Eq. (E3), we have used 


SXA'B'C'D'jq' ,i^n') 
SG c d(k', ie m >) 


= & A' ,c$C' ,dG D' B' {k' - q' liCm'-n') 

3D',c^B*,dG a'C'( k ~\~ q 

(E10) 


where we have used Eq. (44); we have introduced Eq. 
(58) at the final line of Eq. (E3). 

Finally, we can obtain Eqs. |54)-(57) by labeling q 
and correctly as the electron-hole scattering with the 
momentum transfer q and the frequency transfer iQ n . 
For the correct labeling, see Figs. 


Appendix F: Analytic continuations of Eqs. 


In this appendix, we explain the details of the analytic 
continuations of Eqs. (55)—(57) in case that their fre¬ 


quency variables satisfy the inequalities for region 22-11 
or 22-III or 22-IV of TABLE U 

First, we can easily carry out the analytic continuation 
of Eq. ( [55] ) . Namely, since Ime—line' is negative for 
region 22-11 and positive for regions 22-III and 22IV (see 
Table [i]), the MT terms in regions 22-11, 22-III, and 22-IV 
are 


r 22-n;a6c d (^ ^ <?) = ^,0^,0 V^{k k'), (FI) 

r 22-mTa6c d (^ <?) = Ko6 u ,oV™d(k k'), (F2) 


and 


r (l)MT 
^ 22-YV\abcd 


(fc, k'; q) = T 


(1)MT 
22 - 111 ,abed 


(k 7 k’;q), 


(F3) 


respectively. 

Second, we can obtain the AL1 terms for regions 22- 
II, 22-III, and 22-IV after several calculations for the 
analytic continuations by using a similar way for the 
analytic continuations in Sect. II B 1. Replacing the 
sum about bosonic Matsubara frequency by the cor¬ 
responding contour integral 93 and using the analytic 
properties of the single-particle Green’s functions and 
Watd-ABCDiq^^n,] <j 2 ,*^n 2 ), we obtain the AL1 term 
in region 22-11: 


<?' {A} JC 


47T2 


X W aBcA-jCbn(q - q'l ikl n - w'; -q 


-u/;-q',-u/) 


-^EE GcA^k' + q', ie m ’)G BD (k + q ', ie m )W^ BcA . dCbD (q - q', ifl n ; -q', 0) 
q' {A} 


T 

-]vEE G C A(k' + q', ie m ' + iVL n )G BD (k + q', ie m + i^ n )W^ BcA . dCbD (q - q', 0; -q', -ifl n ) 

q' {A} 

tanh UJ ^^ ImG cA( k ' + 9 ')G^,(fc + q'Wa B ^ C bD{q - q'\ -V) 


-^L/‘ 

q' {A} J 

1 

N E- 


+ 


q' {A} 

(Principal integral terms) 


OO 

tanh i ^^ G cA( k ' + q') lmGi sl( k + q'WaBcA^cbDiq - -</) 


(F4) 


where we have used the contour C shown in Fig. H! a ) 


and introduced W^ A ^ CD (q 1 -,q 2 ), 

^a C 6 ( cd } (91, Wl) (q 2 , w 2 : 


- ulLuli 


d r\ • 
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(a) (b) (c) (d) 



FIG. 23: Contours used for the analytic continuations of the AL1 term in (a) region 22-11 and (b) region 22-III or 22-IV and 
of the AL2 term in (c) region 22-11 or 22-III and (d) region 22-IV. 


Furthermore, we have not explicitly written the terms of 
the principal integral since those terms exactly cancel out 
the terms of the principal integral for region 22-III or 22- 
IV; due to this cancellation, such terms are unnecessary 
to calculate the kernel of the CVCs since part of the 


kernel, Jl 


(i) 


22 ;cdCD 


(k, k'\ 0), is proportional to the difference 
betwe en the AL1 terms fo r th e different regions [see Eq. 
(A10)]. In deriving Eq. (F4), we have used a relation 


of the single-particle Green’s function as a result of the 


time-reversal and the even-parity symmetry, G^\k) = 
G^’ik)*; as described in Sect. II C, for more general 
derivation, we should replace ImG^'fl:) by ^[G^\k) — 

G^t\k)]. (The similar remark holding in the other cases 
of the AL term is not described below.) Then, in the 
similar way for region 22-11 by using the contour C shown 
in Fig. mb), we obtain the AL1 term in region 22-III or 
22-IV: 


r (i)ALi n ii \ 1 du' w -t-e ^y(R) / 1 / . , /\w AL ( RA ) 

r 22 -m-abcA k ^ k ;v) = - / vv tanh 

1 q' {A} J-00 


du}' 


2tt 2 T CA 


-jvEE/ tanh 

q' Ml - 

+ (Principal integral terms). 


-^Im+ J)G™{k + q')W:^ cbD {q - g'; ~q') 
—G£l(k' + q')lmG™(k + q')W^% bD {q g'; -g') 


(F6) 


and 


r k q) = rglt^ abcd (k, k 1 - g). (F7) 


Third, carrying out the analytic continuation of the 
AL2 term in the similar way for the AL1 term by using 


the contour C shown in Fig. 23 [c) or|23|(d), we obtain the 
AL2 term for region 22-11 or 22-III or 22-IV. Namely, the 
AL2 terms in regions 22-11, 22-III, and 22-IV are given 
by, respectively, 


t(1)AL2 
22-11 \abcd 


dco 1 


(fc,fc / ;g)=--^^ / tanh 

q, oo 

1 f°° du)' 

— — / / / —— tanh 

N ^ ^ J ^ 2tt 

q, {A} ^-oo 

+ (Principal integral terms), 

p(l)AL2 /, w \ —p(l)AL2 /, w \ 

1 22 -III;abcdl ,C ’ A , gj —1 2 2 -II;abccdV K ’ dh 


2 T 


[—Im G^(k’ - q , )]G^{k + q + q')W^^, cbD {-q'; -q - g') 


i(R) 


AL(RA) 


w + 6 + UJ G^(k' - q')lmG { sl(k + g + q')W^^ cbD (-q'-, ~q - <l') 


2 T 


(F8) 

(F9) 
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and 

p(l)AL2 /, 7 /. \ _ 

1 22-1 V;abcd\ h ’ ^ iV) ~ 


1 (L/ 

- — > > / — tanh ■ 

N J ^ 2it 
r dw' 

— — > > / —— tanh 

N ^ ^ J „ 27r 

q' {A} J 


[-ImG^ik' - g')]G^,(fc + q + “9 ~ ?') 


dA) 


AL(RA) 


2 T 


W + 6 + W G^(fc' - (/c + g + ^)W^ B f d ^ cbD {-ci \ “9 - ?') 


2T 


(Principal integral terms). 


(F10) 


Appendix G: Technical details about numerically 
solving the Bethe-Salpeter equation for the current 


In this appendix, I give several technical remarks about 
the numerical calculations to self-consistently determine 
the current from the Bethe-Salpeter equation. 

Before the remarks about the iterative self-consistent 
procedures of the Bethe-Salpeter equation, I remark on 
how to numerically calculate the frequency integral. We 
calculate the frequency integral by discretizing it with fi¬ 
nite interval such as Ae,- and approximating the upper 
and the lower value of the integral by the finite cutoff val¬ 
ues such as e c and — e c , respectively; those values can be 
appropriately chosen within the numerical accuracy. (In 
the similar way, we numerically calculate the frequency 
integrals of the conductivities.) 

We can self-consistently solve the Bethe-Salpeter equa¬ 
tion for the current including the E CVC, the MT CVC, 
and the AL CVC by iteration in the following procedures. 


First, we determine the MT CVC in the presence 
of the E CVC. Using the input of the current, which 
is A^ 2 . AB (k,€j;0) for the first iteration, we deter¬ 
mine A„ ; 2 ;af>(fe, e?; 0) from Eq. (74). After carrying 
out the Fourier transformations of A„ ; 2 ; cD(fc, ey; 0) and 
ImV5H n (fc — k', e,- — e',) about momentum, we calculate 
AA^ cd ( r ,e,; 0 ) as follows: 


A A mt 

^ A z/;2 :cd 


>>^;0)=EE^f F crfe,4;T) 

C,D e'„ 


x Im Kcir.( r > e j -4) 

xA „ ; 2 jcdC?", e'd 0). (Gl) 


Since F^ T (ej , e',; T) includes the hyperbolic cotangent 
[see Eq. (65)] and its principal integral has the 0/0 struc- 
ture due to hnV^ Cdn (r. 0) = 0, the principal integral can 


be calculated as follows 10 : 


E 


Ae/, 

2n 


coth 


2 T 


ImU 


(R) 


cCdD 


- 4) 


x A „;2 ,cd(t, e',; 0 ) 


= E 


Ac', 

V 

2 tt 


■coth■ 


^ e 


2 T 


J ' lmV cCdD( r ’ e 3 - e 'j') 


x A„ ; 2 ;cc(r, e'j,] 0) 

[( e ^~ + l ) lmV cCdDi. r ^j - 4 


2 71 de' 3 , 


x A„; 2 -,cd{v, 4 ; 0 )] 


(G2) 


After the Fourier transformation of AA^ cd (r, e^; 0) 
about r, we obtain the MT CVC, AA^J. cd (k, ej] 0), and 
then add this to the input of the current. (If we consider 
only the MT CVC, i.e. neglect the AL CVC, we skip 
the following second and third steps, and the sum of the 
input and the MT CVC becomes the output.) 

Second, we turn to the calculation of the AL1 
CVC. Carrying out the Fourier transformation of 

/r>\ ~ 

Im Gcc J (k,e'j,) about k and using A„ ; 2 ;C c(u e j’i 0 ), which 
is the same for the calculation of the MT CVC, we cal¬ 
culate [ImGA] (—fc; e',, w'„) from the equation, 


[ImGA] c (—fc; 4’4") =E e ’ h ' rimG « ) ( r >4+ w i") 

r 

x A I/;2 ; CC (—r, 4 ; 0). (G3) 


After calculating Ncccc(q,Uj) and N^cc* {q,Wj) from 
Eqs. (78) and (79), respectively, and using Eq. (77), we 
determine A' ALr [fc; ej, <*;'■„) from the following equation: 


fC(R)/ 


v c ALl (fc ; £i , 4, ) = Y ^4 Fu hl 4 ,4,4"; T ) 


x F AL1 (ej , 4; T)W c al ( ra ) (fc, - ) 
x [lmGA] c (-fc;e',,4,). (G4) 


In the above equation, we can calculate the principal inte¬ 
gral of the hyperbolic cotangent of F AL 1 (ej, e'.,; T) in the 
similar way for the MT CVC since that principal integral 
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also has the 0/0 structure due to F^ L1 (ej,€j,ujj,r,T) = 0 
[see Eq. ©]■ Carrying out the Fourier transformation 
of X"" about k and using the equation, 

r u/.,, 

3 

x ImG™(r,q+w'„) 
x (G5) 

we obtain the AL1 CVC, A.A^} cc (k, Cj\ 0). At this point, 
we add the AL1 CVC to the sum of the input of the 
current and the MT CVC. 

Third, we calculate the AL2 CVC. Using 
[lmGA] c (fc; e'-,, — w'„) and W^ L ^ RA \k,—ojy,) (which 
have been determined in the above second step), we 
calculate X^ L2 (k\ ej,ujj„) from the following equation: 


X. 


AL2 


(fc; A 4 F u L2 (e 


3 5 C j' 1 “'j 


U;T) 


x i^ L2 (e i; e',; T)TV c AL ( RA )(fc, -wj„) 
x [lmGA] c (fc;4,-w'„). (G6) 

Since F t ^ L2 (ej ,—ej, u )’-„; T) = 0 is satisfied in the above 
equation [see Eq. (69)], we can calculate the principal 
integral of the hyperbolic cotangent of F^ L2 (ej, e',; T) 
in the similar way for the MT CVC. After the Fourier 
transformation of X^ h2 (k; Cj, ojj„) about fc, we obtain 
the AL2 CVC, AA^? ce (fc, e i; 0): 


AA^f cc (fe, e,; 0 )=- 4 ^E ^ E A 4' 

V Cl?'./, 

3 

x V^ L2 (r;ej,<4,)- (G7) 

Adding the AL2 CVC to the sum of the input of the 
current, the MT CVC, and the AL1 CVC, we obtain 
the output of the current for the Bethe-Salpeter equation 
with the E CVC, the MT CVC, and the AL CVC. 

Finally, estimating the sum of the difference between 
the absolute values of the output and the input of the 
current, we judge whether the output can be regarded 
as the solution of the Bethe-Salpeter equation within the 
numerical accuracy: if the difference becomes less than 
10~ 4 , the solution is assumed to be obtained; otherwise, 
after replacing the input by the output, we solve the 
above procedures again. 


Appendix H: Method to determine the dominant 
fluctuations in a multiorbital Hubbard model 

In this appendix, I explain how to determine the domi¬ 
nant fluctuations in a multiorbital Hubbard model among 


the four kinds of fluctuations. 

To discuss the dominant fluctuations in a multiorbital 
Hubbard model, we need to consider charge fluctua¬ 
tions, spin fluctuations, orbital fluctuations, and spin- 
orbital-combined fluctuation;^. This is because the 
Hubbard interaction terms can be expressed in terms 
of their operators^. Then, spin and spin-orbital- 
combined fluctuations are described by Xabcdi^^^n), 
and charge and orbital fluctuations are described by 
Xabcdi.^ i ^nW m - For example, spin fluctuation s are 
characterized by x s (q,ifi n ) = E Q ,6 xlabbi*!’ 
spin-orbital-combined fluctuations are characterized by 
the correlation function between the products of the 
spin and the orbital operator, e.g. E„[xlL:w(< 7 AU„) + 

xL 23 (^)FMi 

Since those four kinds of fluctuations have the differ¬ 
ent dependence on Jp or U, the dominant fluctuations 
can be determined by analyzing Jp and U dependences 
of A^ lax (q) _1 and A[/ lax (< 7 ) _1 . The difference in the J H 
or U dependence arises from the difference in the depen¬ 
dence of the bare four-point vertex function character¬ 
izing the fluctuation. For example, in cas^^ of a t 2 g - 
orbital Hubbard model with J' = Jp and V — U — 2Jp, 
the bare four-point vertex function for charge fluctua¬ 
tions is —U — 4U' + 2Jp = —517 + 10Jp; that for spin 
fluctuations is U + 2Jp; that for orbital fluctuations is 
U' - 2,/p - ,/' = U- 5 Jp or U' - 2 Jp + J' = U - 3 Jp or 
— U + 2U’ — Jp = U — 5Jp; that for orbital-spin-combined 
fluctuations is U' + J' = U — Jp or U' — J' = U — 3 Jp 
or 17 — Jp. Thus, if orbital or spin-orbital-combined fluc¬ 
tuations are dominant, we obtain non-monotonic Jp de¬ 
pendence of A® ax (q) -1 or A laax (q) _1 , respectively, with 
increasing Jp since the bare four-point vertex function 
characterizing those fluctuations changes from repulsive 
to attractive at a critical value of Jp as a function of U; if 
spin fluctuations are dominant, we obtain monotonic Jp 
dependence of A^ lax (q) _1 with increasing Jp due to the 
always repulsive bare four-point vertex function charac¬ 
terizing those. Note that charge fluctuations are always 
suppressed in a realistic set of the Hubbard interaction 
terms. 
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